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1  Executive  Summary 


1.1  Overview  and  Main  Results 

This  report  describes  work  performed  from  25  July  2010  to  19  October  2013  under  AFRL  grant  number 
FA8750-1 1-1-0015  entitled  “Studies  on  Optimized  Assured  Cooperative  Communications”.  We  cover  two 
important  tasks:  (1)  determine  optimal  transmission  power  sequence  for  cooperative  hybrid  automatic-repeat- 
request  (H-ARQ)  relaying  protocol  in  wireless  networks;  and  (2)  design  and  optimize  cooperative  commu¬ 
nication  protocols  by  exploring  all  possible  variations  in  time  and  power  domains.  Some  main  results  are 
summarized  in  the  following: 

•  We  first  addressed  the  fundamental  problem  of  assigning  optimal  transmission  power  sequence  for  co¬ 
operative  hybrid  automatic -repeat-request  (H-ARQ)  relaying  protocol  over  quasi-static  Rayleigh  fading 
channels.  A  closed-form  expression  of  the  average  total  transmission  power  was  obtained  first  for  the 
cooperative  H-ARQ  relaying  protocol,  in  which  the  source  may  use  different  transmission  power  level  in 
different  (re)transmission  rounds.  We  determined  the  optimal  power  sequence  by  minimizing  the  average 
total  transmission  power.  However,  the  closed-form  expression  of  the  average  total  power  consumption 
of  the  cooperative  H-ARQ  protocol  is  complicated  in  general,  so  we  developed  a  simple  approximation  of 
the  average  total  transmission  power  that  is  tight  at  high  SNR.  Based  on  the  asymptotically  tight  approx¬ 
imation,  we  were  able  to  identify  the  sequence  of  power  values  that  minimizes  the  average  total  power 
consumption  of  the  cooperative  H-ARQ  protocol  for  any  given  targeted  outage  probability.  In  particular, 
we  derived  a  set  of  equations  that  describe  the  optimal  power  level  in  each  (re)transmission  and  enable  its 
recursive  calculation  with  fixed  searching  complexity.  When  the  maximum  number  of  (re)transmissions 
allowed  in  the  protocol  is  L  =  2,  we  have  a  closed-form  result  for  the  optimal  transmission  power  se¬ 
quence.  The  optimal  power  assignment  solution  reveals  that  conventional  equal  power  assignment  is 
not  optimal  in  general.  For  example,  for  targeted  outage  probability  of  I  O  '4  with  a  maximum  of  two 
(re)transmissions,  the  average  total  transmission  power  with  the  optimal  power  assignment  is  3  dB  lower 
than  the  equal  power  assignment.  The  difference  in  average  total  power  cost  grows  further  when  the  num¬ 
ber  of  allowable  (re)transmissions  increases  (for  example,  6  dB  gain  with  a  cap  of  4  (re)transmissions)  or 
the  targeted  outage  probability  decreases  (11  dB  gain  with  outage  probability  10-5  and  (re)transmission 
capped  at  4). 

•  Cooperative  communication  has  emerged  as  a  new  wireless  network  communication  concept,  in  which 
parameter  optimization  such  as  power  budget  and  time  allocation  plays  an  important  role  in  cooperative 
relaying  protocol  designs.  While  most  existing  works  on  cooperative  relaying  protocol  designs  consid¬ 
ered  equal-time  allocation  scenario,  i.e.  equal  time  duration  is  assigned  to  each  source  and  each  relay,  in 
this  project  we  designed  and  optimized  cooperative  communication  protocols  by  exploring  all  possible 
variations  in  time  and  power  domains.  First,  we  considered  an  ideal  cooperative  relaying  protocol  where 
the  system  can  use  arbitrary  re-encoding  methods  at  the  relay  and  adjust  time  allocation  arbitrarily.  We 
obtained  an  optimum  strategy  of  power  and  time  allocations  to  minimize  the  outage  probability  of  the 
ideal  cooperative  protocol.  Specifically,  for  any  given  time  allocation,  we  were  able  to  determine  the 
corresponding  optimum  power  allocation  analytically  with  a  closed-form  expression.  We  also  showed 
that  in  order  to  minimize  the  outage  probability  of  the  protocol,  one  should  always  allocate  more  energy 
and  time  to  the  source  than  the  relay.  Second,  with  more  realistic  consideration,  we  designed  a  practi¬ 
cal  cooperative  relaying  protocol  based  on  linear  mapping,  i.e.  using  linear  mapping  as  the  re-encoding 
method  at  the  relay  and  considering  integer  time  slots  in  the  two  phases.  The  theoretical  results  from 
the  ideal  cooperative  protocol  served  as  guideline  and  benchmark  in  the  practical  cooperative  protocol 
design.  We  also  developed  an  optimum  linear  mapping  to  minimize  the  outage  probability  of  the  linear- 
mapping  based  cooperative  protocol.  Extensive  numerical  and  simulation  studies  illustrate  our  theoretical 
developments  and  show  that  the  performance  of  the  proposed  cooperative  relaying  protocol  based  on  the 
optimum  linear  mapping  is  close  to  the  performance  benchmark  of  the  ideal  cooperative  protocol. 
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The  introduction  to  the  project  and  an  overview  of  each  research  task  is  given  in  Section  2.  Then  in  Section 
3  we  present  specifically  models,  assumptions,  methods,  and  main  results  for  each  task.  Section  4  contains  the 
conclusions  and  bibliography  is  included  at  the  end. 
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2  Introduction  to  the  Project 

2.1  Overview  of  Cooperative  ARQ  Relaying  Protocols  in  Wireless  Networks 

Cooperative  wireless  networks  can  substantially  increase  network  reliability  as  each  user’s  information  may 
be  jointly  delivered  to  its  destination  with  the  assistance  of  cooperative  users/nodes  in  the  networks  [1] — [8]. 
On  the  other  hand,  automatic-repeat-request  (ARQ)  protocols  have  been  long-time  used  to  enable  reliable  data 
packet  transmissions  in  data  link  control  [9]— [12],  in  which  a  receiver  requests  retransmission  when  a  packet 
is  not  correctly  received.  In  basic  ARQ  protocols,  a  receiver  decodes  an  information  packet  based  only  on 
the  received  signal  in  each  (rc)transmission  round  [9,  10],  while  in  more  advanced  ARQ  protocols,  a  receiver 
may  decode  an  information  packet  by  combining  received  signals  from  all  previous  (re)transmission  rounds, 
resulting  in  so-called  hybrid  ARQ  (H-ARQ)  protocols  [11,  12]. 

It  is  a  natural  idea  to  exploit  H-ARQ  protocols  in  conjunction  with  the  cooperative  communication  concept  to 
jointly  enhance  link  connectivity  and  network  reliability.  When  a  source  sends  a  wireless  signal  to  an  intended 
destination,  nearby  users  may  also  receive  the  signal.  Thus,  if  the  destination  requests  retransmission,  the  nearby 
users  may  also  assist  forwarding  the  signal  alongside  the  source’s  retransmission.  The  destination  can  combine 
all  received  signals  from  source  and  relays  and  jointly  decode  the  signals  to  improve  detection  performance. 
Some  recent  works  have  studied  H-ARQ  protocols  in  the  context  of  cooperative  relay  networks  [13]— [16].  In 
particular,  [13]  was  among  the  first  such  studies  to  present  a  general  framework  of  cooperative  ARQ  relay  net¬ 
works.  It  was  shown  that  cooperative  ARQ  relay  networks  have  great  advantages  in  terms  of  throughput,  delay, 
and  energy  consumption  compared  to  conventional  multihop  ARQ  networks  in  which  point-to-point  ARQ  links 
are  concatenated  to  form  network  routes.  In  [14,  15],  information-theoretic  analysis  was  developed  and  upper 
bounds  for  the  diversity  order  of  a  decode  -  a  n  d-  fo  rw  aid  (DF)  cooperative  ARQ  relay  scheme  were  characterized 
for  both  slow  and  fast  fading  channels  as  a  means  to  study  the  diversity-multiplexing-delay  tradeoff.  Recently, 
a  closed-form  asymptotically  tight  approximation  of  the  outage  probability  for  the  DF  cooperative  ARQ  re¬ 
lay  scheme  was  developed  under  fast  fading  conditions  [16].  It  was  shown  that  the  cooperative  ARQ  scheme 
achieves  diversity  order  equal  to  2L  —  1  while  the  diversity  order  of  the  direct  ARQ  transmission  scheme  is 
only  L,  where  L  is  the  maximum  number  of  (re)transmission  rounds. 

In  wireless  links  formed  by  wireless  devices  with  limited  power  resources,  power  efficiency  is  a  key  research 
matter  in  the  optimization  of  ARQ  retransmission  protocols  [17] — [21].  In  both  [17]  and  [18],  the  power  ef¬ 
ficiency  of  ARQ  protocols  was  examined  under  the  assumption  of  the  same  transmission  power  level  in  each 
(re)transmission  round.  In  [19],  by  assuming  that  partial  channel  state  information  (CSI)  is  available,  optimal 
transmission  power  in  each  (re)transmission  round  was  determined  for  an  H-ARQ  protocol  by  a  linear  pro¬ 
gramming  method  that  selects  a  power  value  from  a  set  of  discrete  power  levels.  In  [20],  without  assuming 
CSI  available  at  the  transmitter  side,  an  optimal  power  transmission  strategy  was  identified  for  a  basic  ARQ 
protocol  where  the  receiver  decodes  based  only  on  the  received  signal  in  each  (re)transmission  round.  It  was 
assumed  that  the  channel  changes  independently  in  each  round.  In  [21],  by  a  recursive  calculation,  an  optimal 
power  assignment  sequence  for  an  H-ARQ  protocol  was  determined  in  quasi-static  Rayleigh  fading,  in  which 
the  channel  does  not  change  during  (re)transmissions  of  the  same  information  packet. 

Note  that,  while  the  average  total  power  consumption  needed  in  the  delivery  of  each  information  packet  has 
been  well  understood  for  the  conventional  non-cooperative  H-ARQ  protocols  [19] — [21],  the  study  of  coopera¬ 
tive  H-ARQ  counterpart  has  been  proven  very  challenging  with  no  available  results  until  this  present  work.  In 
this  project,  we  try  to  determine  the  optimal  power  sequence  assignment  that  minimizes  the  average  total  power 
consumption  for  the  cooperative  H-ARQ  relaying  protocol  in  which  each  relay  forwards  Alamouti-based  re¬ 
transmission  signals.  We  consider  a  quasi-static  Rayleigh  fading  environment  by  assuming  that  the  channels  do 
not  change  during  (re)transmissions  of  the  same  information  packet  and  they  may  change  independently  when 
the  protocol  transmits  a  new  information  packet.  We  develop  a  new  analytical  approach  and  obtain  a  closed- 
form  expression  of  the  average  total  transmission  power.  The  closed-form  result  is  valid  for  any  maximum 
number  of  (re)transmission  rounds  L  allowed  in  the  protocol  in  which  the  source  may  use  varying  transmission 
power  level  per  round.  The  closed-form  expression  may  serve,  therefore,  as  the  basis  to  optimize  power  alloca¬ 
tion  for  the  cooperative  H-ARQ  relaying  protocol.  However,  the  closed-form  expression  is  rather  complicated 
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for  large  L,  so  we  develop  a  simple  approximation  of  the  average  total  transmission  power  which  is  tight  at  high 
SNR.  Based  on  the  asymptotically  tight  approximation,  we  arc  able  to  identify  the  sequence  of  power  values 
that  minimizes  the  average  total  power  consumption  of  the  cooperative  H-ARQ  relaying  protocol  for  any  given 
tai'gctcd  outage  probability.  In  particular,  we  derive  a  set  of  equations  that  describe  the  optimal  power  level  in 
each  (re)transmission  and  enable  its  recursive  calculation  with  fixed  searching  complexity.  When  the  maximum 
number  of  (re)transmissions  allowed  in  the  protocol  is  L  =  2,  we  have  a  closed-form  result  for  the  optimal 
transmission  power  sequence.  The  optimal  power  assignment  solution  reveals  that  conventional  equal  power 
assignment  scheme  which  uses  the  same  transmission  power  in  all  (re)transmission  rounds  is  not  optimal.  As 
an  example,  for  a  targeted  outage  probability  of  10“  4  and  the  maximum  number  of  (re)transmissions  L  =  3, 
the  average  total  transmission  power  based  on  the  optimal  power  assignment  is  6  dB  less  than  that  of  using 
the  common  equal  power  scheme.  We  also  observe  that  the  larger  the  maximum  number  of  (re)transmissions 
allowed  in  the  cooperative  H-ARQ  protocol  or  the  lower  the  required  outage  probabilities,  the  more  power 
savings  the  optimal  power  assignment  strategy  offers.  Extensive  simulation  and  numerical  results  arc  provided 
to  illustrate  and  validate  our  theoretical  development. 

2.2  Overview  of  Cooperative  Communication  Protocol  Designs  Based  on  Optimum  Power  and 
Time  Allocation 

Cooperative  communication  has  received  significant  attention  recently  as  an  emerging  communication  concept 
for  wireless  networks  [24,  25].  Due  to  the  broadcast  nature  of  wireless  transmissions,  cooperative  communica¬ 
tions  enables  neighboring  network  nodes  to  share  resources  and  cooperate  to  send  information  to  an  intended 
node.  Distributed  transmissions  from  source  and  relay  nodes  provide  spatial  diversity  (as  well  as  multiplex¬ 
ing  gain  in  some  designs)  for  information  detection  at  a  destination  node.  Cooperative  communications  can 
significantly  improve  system  performance  and  robustness  of  wireless  networks  especially  in  severe  fading  en¬ 
vironment.  Cooperative  relaying  techniques  have  been  considered  in  some  latest  communication  standards,  for 
example,  in  IEEE  802. 16j  WiMAX  standard  [26]  and  3GPP’s  Long  Term  Evolution  (LTE) -Advanced  standard 
[27], 

The  idea  of  cooperative  communications  can  be  traced  back  to  1970’s  [28]  and  information-theoretic  studies 
have  been  extensively  carried  out  since  then  (see  [29,  30,  31,  32,  33]  and  the  references  therein).  In  recent 
years,  besides  information-theoretic  studies,  many  efforts  have  been  shifted  to  design  practical  cooperative 
communication  protocols  for  wireless  networks  with  specific  system  constraints  and  quality-of-service  (QoS) 
requirements.  Various  cooperative  relaying  protocols  such  as  dec  ode  -  a  nd-  fo  rw  aid  (DF)  relaying  protocol  and 
amplify-and-forward  (AF)  relaying  protocol  were  proposed  for  wireless  networks  (see,  for  examples,  [34,  35, 
36,  37])  and  substantial  performance  gains  of  such  relaying  protocols  have  been  demonstrated  compared  to 
conventional  non-cooperative  transmission  approach.  Cooperative  relaying  protocols  have  been  generalized  to 
multi-relay  networks  with  either  parallel  or  sequential  multiple  relays  to  further  improve  network  performance 
with  price  of  higher  complexity  and  more  overhead  [34,  38].  Cooperative  communication  protocols  have  also 
been  integrated  into  conventional  QoS  control  mechanisms  such  as  automatic -repeat-request  (ARQ)  protocols 
to  enhance  reliability  and  robustness  of  wireless  networks  [39,  40].  Cooperative  relaying  techniques  have  been 
applied  to  multimedia  wireless  broadcast  and  multicast  applications  with  substantial  data  rate  increase  and 
power  saving  [41,  42].  More  thorough  discussions  on  basic  theories,  protocols,  and  applications  for  cooperative 
communications  can  be  found  in  [24,  25]. 

Resource  allocation  in  cooperative  communication  protocols  such  as  power  budget  and  time  allocation  to 
source  and  relays  plays  an  important  role  in  the  overall  performance  of  the  protocols.  Note  that  most  existing 
works  on  cooperative  relaying  protocol  designs  considered  equal-time  allocation  scenario,  i.e.  equal  time  du¬ 
ration  is  assigned  to  each  source  and  each  relay  [34,  35,  36,  37,  38,  39,  40,  41,  42,  43,  44,  45].  For  example, 
with  equal-time  allocation,  in  [43]  source  and  relay  power  allocation  was  optimized  such  that  the  outage  prob¬ 
ability  of  the  cooperative  relaying  protocol  is  minimized.  In  [44],  also  with  equal-time  allocation,  optimum 
power  allocation  was  determined  in  which  symbol  error  rate  (SER)  performances  of  both  DF  and  AF  relaying 
protocols  were  optimized  respectively.  However,  there  are  rather  limited  studies  on  non-equal  time  allocation 
scenario  and  further  on  joint  optimization  of  power  allocation  and  time  allocation  in  cooperative  relaying  pro- 
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tocol  designs.  In  [46],  by  assuming  that  instantaneous  channel  state  information  (CSI)  of  source-destination, 
source-relay  and  relay-destination  links  arc  available,  joint  power  and  time  allocation  were  optimized  in  order 
to  minimize  the  outage  probability  of  the  cooperative  communication  system.  In  [47],  it  was  shown  that  the 
problem  of  minimizing  the  outage  probability  with  respect  to  joint  power  and  time  allocation  at  source  and 
relays  is  a  multi-variable  convex  optimization  problem  for  nonorthogonal  cooperative  communication  systems 
where  the  source  and  the  relay  arc  allowed  to  transmit  simultaneously.  However,  numerical  method  was  consid¬ 
ered  in  [47]  to  solve  the  convex  optimization  problem  without  analytical  solutions.  It  is  hard  to  obtain  insight 
understanding  of  the  cooperative  relaying  protocols  based  on  the  numerical  results.  We  note  that  in  both  [46] 
and  [47],  independent  codebooks  were  used  for  source  and  relays  in  the  cooperative  relaying  protocol  designs. 

In  this  project,  we  intend  to  design  and  optimize  cooperative  communication  protocols  by  exploring  all 
possible  variations  in  time  and  power  domains.  We  assume  that  the  channels  are  known  at  the  receiver  side, 
but  not  at  the  transmitter  side.  We  consider  an  orthogonal  cooperative  communication  scenario,  in  which  the 
source  transmits  signals  in  Phase  I  and  the  relay  decodes  the  signals  and  forwards  them  to  the  destination  in 
Phase  II.  First,  we  analyze  an  ideal  cooperative  communication  protocol  where  the  system  can  use  arbitrary 
re-encoding  methods  at  the  relay  and  adjust  time  allocation  arbitrarily  between  Phases  I  and  II.  We  are  able 
to  obtain  an  optimum  strategy  of  power  and  time  allocations  to  minimize  the  outage  probability  of  the  ideal 
cooperative  protocol.  Specifically,  for  any  given  time  allocation,  we  arc  able  to  determine  the  corresponding 
optimum  power  allocation  at  the  source  and  relay  analytically  with  a  closed-form  expression.  We  also  show 
theoretically  that  in  order  to  minimize  the  outage  probability  of  the  protocol,  one  should  always  allocate  more 
energy  and  time  to  Phase  I  than  that  to  Phase  II  in  the  protocol. 

Note  that  in  the  ideal  cooperative  protocol  in  which  there  is  no  constraint  on  the  re-encoding  methods  and 
time  allocation,  it  may  not  be  easy  or  feasible  to  implement  it  in  practical  systems.  So,  next  we  would  like 
to  propose  a  practical  cooperative  communication  protocol  design  based  on  linear  mapping.  The  theoretical 
results  from  the  ideal  cooperative  protocol  will  serve  as  guideline  and  benchmark  in  the  practical  cooperative 
protocol  design.  More  specifically,  we  design  a  cooperative  communication  protocol  by  considering  linear 
mapping  as  the  re-encoding  method  at  the  relay,  where  the  protocol  uses  integer  time  slots  in  Phases  I  and 
II.  It  is  much  easier  to  implement  the  linear  mapping  forwarding  method  with  the  time  allocation  of  integer 
time  slots  in  the  two  phases  in  practical  systems.  We  also  develop  an  optimum  linear  mapping  to  minimize 
the  outage  probability  of  the  linear-mapping  based  cooperative  protocol.  Interestingly,  simulation  results  show 
that  the  performance  of  the  proposed  cooperative  protocol  based  on  the  optimum  linear  mapping  is  close  to  the 
performance  benchmark  of  the  ideal  cooperative  protocol. 


3  Models,  Assumptions,  Methods  and  Procedures 

3.1  Cooperative  Hybrid-ARQ  Relaying  in  Wireless  Networks:  Average  Power  Consumption 
and  Optimal  Power  Assignment 

In  this  subsection,  we  determine  optimal  transmission  power  sequence  for  cooperative  hybrid  automatic-repeat- 
request  (H-ARQ)  relaying  protocol  in  wireless  networks.  First,  we  describe  briefly  the  cooperative  H-ARQ 
relaying  protocol.  Second,  we  derive  the  average  total  transmission  power  expended  in  the  delivery  of  each 
information  packet.  Third,  we  develop  an  approximation  of  the  average  total  transmission  power  which  is  tight 
at  high  SNR.  Fourth,  we  determine  an  optimal  power  assignment  strategy  for  the  cooperative  H-ARQ  relaying 
protocol  based  on  the  approximation  of  the  average  total  transmission  power  developed  in  the  previous  section. 
Finally,  numerical  and  simulation  studies  arc  carried  out  to  compare  the  performance  of  the  optimal  and  equal 
power  assignment  strategies. 

3.1.1  System  Model 

For  simplicity  in  presentation,  we  consider  a  cooperative  H-ARQ  relaying  model  with  one  source,  one  relay 
and  one  destination,  as  illustrated  in  Fig.  1.  We  assume  that  L  is  the  maximum  number  of  (re)transmission 
rounds  allowed  in  the  protocol.  The  cooperative  H-ARQ  relaying  scheme  operates  as  follows.  First,  the  source 
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Figure  1 :  Illustration  of  a  cooperative  H-ARQ  relaying  protocol  with  source  (re)transmission  power  Ps,m  ( 1  < 
m  <  L)  and  relay  power  Pr>n  (k  +  1  <  n  <  L)  when  the  relay  decodes  correctly  at  the  fc-th  round  and  stalls 
cooperating  at  round  (k  +  1). 


broadcasts  an  information  packet  to  the  destination  and  the  relay.  The  destination  sends  a  single  bit  of  acknowl¬ 
edgement  (ACK)  or  negative-acknowledgement  (NACK)  indicating  success  or  failure  of  receiving  the  packet, 
respectively,  to  both  the  source  and  the  relay.  The  ACK/NACK  feedback  is  assumed  to  be  detected  error-free 
at  the  source  and  the  relay.  If  ACK  is  received  by  the  source  or  retransmission  reaches  the  maximum  number 
of  rounds  L,  the  source  begins  transmission  of  a  new  information  packet.  If  NACK  is  received  by  the  source 
and  the  maximum  number  of  rounds  L  is  not  reached,  the  source  retransmits  the  packet  at  a  potentially  differ¬ 
ent  transmission  power.  If  the  relay  decodes  successfully  ahead  of  the  destination,  the  relay  stalls  cooperating 
with  the  source  by  forwarding  the  packet  to  the  destination  by  using  a  space-time  transmission  scheme  [14], 
for  example  the  Alamouti  scheme  [23].  The  destination  combines  the  signals  from  the  source  and  the  signals 
from  the  relay  and  jointly  decodes  the  information  packet.  If  the  destination  still  cannot  decode  an  information 
packet  after  L  (re)transmission  rounds,  an  outage  event  is  declared  which  means  that  the  signal-to-noise  (SNR) 
of  the  combined  received  signals  is  below  a  required  SNR.  The  probability  of  an  outage  event  is  referred  to  as 
the  outage  probability. 

The  cooperative  H-ARQ  relaying  scheme  can  be  modeled  as  follows.  The  received  signal  yr>m  at  the  relay  at 
the  rn-th  (1  <  m  <  L)  H-ARQ  (re)transmission  round  can  be  modeled  as 

Ur,m .  =  \J  Ps,mhsrXs  +  hr, mi  (1) 

where  Ps,m  is  the  source  transmitted  power  at  the  m-th  H-ARQ  (re)transmission  round,  hsr  is  the  coefficient 
of  the  source-relay  channel,  xs  is  the  transmitted  information  packet  from  the  source,  and  r]r,m  is  additive 
noise.  If  the  relay  is  not  involved  in  forwarding,  the  received  signal  y,i.rn  at  the  destination  at  the  m-th  H-ARQ 
(re)transmission  round  is 

Vd,m  —  \/ Ps,mhsd%s  T  hd,mi  (2) 

where  hsd  is  the  source-destination  channel  coefficient. 

If  the  relay  decodes  the  packet  from  the  source  successfully,  it  helps  in  forwarding  it  to  the  destination  using 
the  Alamouti  scheme.  It  is  assumed  that  the  relay  knows  the  codeword  of  the  packet.  Specifically,  if  the  packet  is 
partitioned  into  two  parts  as  xs  =  [.xs.j  xsp\,  then  the  relay  forwards  a  corresponding  vector  xr  =  [— x*  2  x*  ,  ] . 
The  received  signal  y(]  .m  at  the  destination  at  the  m-th  H-ARQ  (re)transmission  round  with  relay  forwarding 
can  be  written  as 

Vd,m  —  yj Ps,mhsd^ s  T  \/ Prhrd%r  T  hd,mi 
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(3) 


where  Pr  is  the  relay  transmitted  power  which  is  assumed  to  be  fixed  over  all  retransmission  rounds1  and  hrd/ 
is  the  channel  coefficient  from  the  relay  to  the  destination.  At  the  destination,  the  packet  xs  can  be  recovered 
based  on  the  orthogonal  structure  of  the  Alamouti  code.  The  destination  combines  the  received  signals  from  all 
(re)transmission  rounds  and  jointly  decodes  the  information  packet.  We  consider  quasi-static  Rayleigh  fading 
channels,  i.e.,  the  channel  coefficients  arc  assumed  to  be  fixed  during  (re)transmissions  of  a  packet  and  may 
change  independently  when  a  new  information  packet  is  transmitted.  The  channel  state  information  is  assumed 
to  be  known  at  the  receiver  side  and  unknown  at  the  transmitter  side.  The  channel  coefficients  hS([,  hsr  and  hrd/ 
arc  modeled  as  independent,  zero-mean  complex  Gaussian  random  variables  with  variances  1  / Asd,  1  / Xsr  and 
1  / \rc],  respectively.  The  noise  variances  r/r.m  and  r/,/  m  are  modeled  as  zero-mean  complex  Gaussian  random 
variables  with  variance  A/q. 


3.1.2  Average  Total  Transmission  Power  of  the  Protocol 

In  this  subsection,  we  first  study  the  overall  SNR  at  the  relay  and  at  the  destination  resulting  from  the  cooperative 
H-ARQ  retransmissions.  Then,  we  calculate  the  probability  of  the  event  that  the  cooperative  H-ARQ  relaying 
protocol  stops  successfully  at  the  f-th  (1  <  l  <  L)  round.  Finally,  we  arc  able  to  derive  the  average  total 
transmission  power  of  the  cooperative  H-ARQ  protocol. 

The  relay  is  to  combine  the  received  signals  from  all  previous  (re)transmission  rounds  from  the  source  based 
on  maximal  ratio  combining  (MRC)  technique  [22]  and  jointly  decode  the  information  packet.  Then,  the  SNR 
of  the  MRC  output  at  the  relay  at  the  l- th  (1  <  l  <  L)  round  is 


7 'r,l  = 


Em=  1  Ps,m\hs 


The  destination  combines  the  received  signals  from  all  previous  rounds  from  the  source  and  the  relay  by  using 
the  MRC  technique  and  then  jointly  decodes  the  information  packet.  To  determine  the  overall  SNR  of  the 
combined  signal  at  the  destination,  we  need  the  following  lemma. 

Lemma  1  In  the  cooperative  H-ARQ  relaying  scheme,  if  the  relay  decodes  successfully  a  packet  from  the 
source  at  the  k-th  (1  <  k  <  L)  (re)transmission  round  and  starts  forwarding  at  round  k  +  1,  then  the  overall 
SNR  at  the  destination  at  the  l-th  (k  <1  <  L)  round  is  given  by 


ld,l,k  = 


Em=i  Ps,m\hsd\2  +  (l~  k)Pr\hrd p 


Proof  :  First,  let  us  focus  on  a  single  transmission  round.  Since  the  packet  from  the  source  can  be  writ¬ 
ten  as  xs  =  [x.5,i  xsp\  and  the  relay  forwards  xr  =  [  -x*  2  x*  ,  ] ,  so  the  corresponding  received  signals 
[yd,m, i  Dd,m, 2]T  at  the  destination  at  the  m-th  (1  <  m  <  L )  round  are 


Ud,m,  1  \/  Ps,m%s,  1  \J  Pr,mxs,  2  h-sd  Vd,m,  1 

Dd,m, 2  pj Ps,m%s,2  \/ Pr,m%s  \  _ 


where 


p  — 

1  r,m  — 


0,  if  m  <  k; 
Pr,  if  m  >  k. 


(Recall  that  the  relay  is  assumed  to  decode  correctly  at  the  fe-th  round  and  starts  forwarding  at  the  (k  +  l)-th 
round  with  fixed  transmission  power.) 

'For  simplicity  in  the  treatment  herein  as  well  as  practical  system  implementation  purposes,  it  is  assumed  that  relay  helpers  have 
fixed,  non  varying,  power  level  across  retransmissions. 
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The  received  signal  vector  can  be  rewritten  as 


Vd,m,  1 

\J  P 5  ,mhsd 

\J  Pr,mhrd 

Xs,l 

Vd,m,l 

}ld,m,2_ 

_y/Pr,mh*d 

\/Ps,mh*d  _ 

.<  2_ 

1 

Rd,m,2_ 

A Y 
— 1  m 


=Hn 


=N„ 


Thus,  the  SNR  at  the  destination  at  the  m-tli  round  is 


SNRm 


(6) 


where  \\Hnix\\2F  is  the  Frobenius  norm  (or  total  power)  of  Hmx.  Note  that 

\\Hmx\\p  =  tr  (xxnH^Hm)  =  (|xS!i|2  +  |x5j2|2)  ( Ps,m\hsd\ 2  +  Pr,m\hrd\ 2)  • 

Let  us  assume  that  the  length  and  (normalized)  power  of  the  packet  xs  is  J .  Then,  we  have  \\Hrnx\\2F 
J  (Ps,m\hsd\2  +  Pr,m\hrd\2)  and  ||Wm|||,  =  JMq.  Thus,  the  SNR  at  the  destination  at  the  m-th  round  is 


SNRm 


Ps,m\hsd\  +  Pr,m\hrd\ 

A/q 


(V) 


where  Pr,m  =  0  if  m  <  k  and  Pr,m  =  Pr  if  m  >  k.  Therefore,  the  overall  SNR  at  the  destination  by  combining 
all  received  signals  from  the  first  l  rounds  is 


i 

Id, l,k  =  ^2 
m=  1 


SNRm  = 


k 

E 

m= 1 


P s,m|^sd| 

A/n 


l 

+  E 

m=k+ 1 


A/n 


(8) 


which  leads  to  the  result  in  the  Lemma. 


□ 


3.1.3  Probability  of  the  Protocol  Successfully  Completing  at  the  /-th  (1  <  l  <  L)  Round 

When  the  destination  successfully  decodes  an  information  packet,  the  cooperative  H-ARQ  protocol  stops  and 
that  may  happen  at  any  round.  In  this  subsection,  we  derive  the  probability  of  the  event  that  the  protocol  stops 
at  the  /- th  (1  <  l  <  L)  round.  This  probability  is  a  key  result  in  determining  the  average  total  transmission 
power. 

Let  {Tr  =  A:}  denote  the  event  that  the  relay  decodes  successfully  at  the  A:-th  round  and  starts  forwarding  at 
round  {k  +  1),  for  any  k  =  1, 2,  •  ■  ■  ,  l  —  1.  Especially,  let  {Tr  =  1}  denote  the  event  that  the  relay  decodes 
unsuccessfully  in  the  first  l  —  1  rounds  (in  this  case,  the  relay  has  no  participation  in  message  forwarding  when 
the  protocol  finishes  at  the  l- th  round).  Furthermore,  let  1  denote  the  probability  that  the  H-ARQ  protocol 
stops  at  the  first  transmission  round  (l  =  1),  which  means  that  the  destination  decodes  successfully  at  round 
one.  When  2  <  l  <  L  —  1,  let  qij.  denote  the  conditional  probability  that  the  H-ARQ  protocol  stops  at  the  l- th 
round  given  that  the  event  {Tr  =  A;}  occurred.  In  other  words,  qu-  is  the  probability  of  H-ARQ  stopping  at  the  l- 
th  round  under  the  condition  that  the  relay  started  forwarding  at  the  (A:  +  l)-th  round  for  any  k  =  1,  2,  —  1. 

Especially,  let  qi  i  (2  <  l  <  L  —  1)  denote  the  probability  that  the  protocol  stops  successfully  at  the  /-th  round 
before  the  relay  can  help.  Moreover,  when  /  =  L,  let  q^j-  denote  the  conditional  probability  that  the  protocol 
stops  at  the  L-th  round  no  matter  whether  decoding  at  the  L-th  round  is  successful  or  not  under  the  condition 
that  the  relay  started  forwarding  at  the  ( k  +  l)-th  round  for  any  Ar  =  1,2,---  ,  L  —  1.  Especially,  let  denote 
the  probability  that  the  protocol  stops  at  the  last  round  regardless  of  successful  decoding  or  not  (in  this  case 
again  the  relay  has  no  chance  to  help). 

With  the  above  notation,  the  average  total  transmission  power  of  the  cooperative  H-ARQ  relaying  protocol 
can  be  given  by 

L  l 

P(L)  =  Ps,lPstop'1  +  EE  qi,k Pr  [Tr  =  k]  Riyk, 

1=2  k= t 
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(9) 


where  R/  j.  is  the  transmission  power  that  the  source  and  the  relay  spend  totally  up  to  the  l- th  round,  which  is 
given  by 

Rl,k  =  +  (l  ~  k)Pr,  (10) 


in  which 


(ii) 


m=l 


First,  we  calculate  the  probability  Pr  [Tr  =  k]  in  (9).  For  any  given  targeted  SNR  70,  the  probability  that  the 
relay  decodes  the  packet  successfully  at  the  first  round  (Tr  =  1)  is 


Pr  \Tr  =  1]  =  Pr 


Ps,i\hs 

A/n 


>  70 


=  e 


(12) 


where  70  =  70A/0.  For  any  Tr  =  k,  k  =  2,  3, I  —  1,  the  probability  that  the  relay  decodes  successfully  at 
the  A'-tli  round,  which  means  that  the  overall  SNR  is  below  the  targeted  SNR  70  until  the  [k  —  l)-th  round  but 
above  70  at  the  subsequent  k- th  round,  is  calculated  as 

Pr  [Tr  =  k]  =  Pr  [7^-1  <  70, 7 r,k  >  7o] 

P.s.k  |  hSr  | 


=  Pr 


7o 


A/n 


^  7r,fc— 1  70 


Pr  [7r,fc_i  <  70]  -  Pr  [7 r,k  <  7o] 

*sr7Q  *sr% 


=  e  —  6  ^k—i 


(13) 


When  Tr  =  l,  we  have 


Pr  [Tr  =  l]=  Pr  [lr,l-i  <  7o]  =  1  -  e 


^sr7Q 

1 


(14) 


Next,  we  calculate  the  conditional  probability  qiq.  in  (9).  When  l  =  1,  the  probability  that  the  protocol  stops 
at  the  first  round  is 


P 


stop,  1  _  pr 


Ps,i\hsd\" 

A/n 


>  7o 


sd^O 


=  e 


(15) 


When  2  <  l  <  L  —  1,  we  consider  two  scenarios:  (i)  k  =  1,  2, —  1,  and  (ii)  k  =  l.  For  any  k  =  1,  2, —  1, 
the  conditional  probability  qi 7.  is  given  by 


where 


qi,k  =  Pr  hd,l-i,k  <  7o,7 d,i,k  >  7o]  =  Pr  [a  <  \hsd\2  <  b ] 


70  -  (/  -  k)Pr\hrd\2  ^0-(l-k-l)Pr\hrd\2 

a  = - — - ,  b  =  - 


(16) 


0/ 


1- 1 


We  observe  that  b  should  not  be  negative,  i.e., 


b  =  70  ~  (l  ~  k  -  l)Pr\hrd\2  >  Q 


n 


1-1 


which  implies  that  \hrd\2  <  jjzrfrzyypr-  When  a  <  0,  it  means  \hrd\2  >  ■  When  a  >  0,  it  means 

\hrd\2  <  pjfcjp  •  Therefore,  the  conditional  probability  (16)  can  be  calculated  as  follows 


Ql,k  = 


I'd 


70 


( l-k)Pr  I  /  (l-k-l)Pr 

/  +  / 


\_J\hrd\2=0  Ahsd\2=a  J \hrd\2=  (l_fc°)pr  4|fc3d|2=0 

Ai  +  A2  —  As, 


Nd  I  hrd,  1 2 


d\hsd\2d\Kdf 


(17) 
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where 


Ai  = 


f(l  k)Pr  e  A^aArde  Xrd\hrd\2 d\hrd\ 

Jo 


^sd{l  k)P r  ^ 


rd 


Xrd% 


e  «-k)pr  -e  ni  , 


(18) 


_70_ 


Ao  = 


1  ’  r  \rde  xrd\Kd? t\\h  A2 


f  70 
(l-k)Pr 

rd^Q 


d|  ^r<i| 


rd^O 


g  (l-k)Pr  —  g  (Z  — /c  — l)Pr 

^0 


As  = 


{l  k~1)Pr  p-KdbXrde-Kd\Kd\2A\h  ,|2 


/o 


d|Zlrdh 


Xrd 


fil-1 


r*d 


Ard^O 


_ A sd'VO 

g  (l  —  k  —  l)Pr  —  g  ^Z-l 


(19) 


(20) 


In  case  of  k  =  l,  the  conditional  probability  qu  is  given  by 

qi,l  =  Pr  <  lo,ld,i,i  >  7o] 

D  |7o  .  I,  ,2  ,  To 

=  T  “  1  ",l  <77 

sd70  _•  xsd7 0 


=  e  S11  —  e  ‘-1 


(21) 


When  l  =  L,  to  calculate  the  conditional  probability  qjjd;,  we  also  consider  two  scenarios:  (i)  k  =  1,  2, L— 
1,  and  (ii)  k  =  L.  For  any  fc  =  1,2, L  —  1,  the  conditional  probability  q^q.  is  given  by 


where 


QL,k  =  Pr  [7d,L-i,fc  <  To]  =  Pr  [| hsd\2  <  c] 
To  -  (L-k  -  l)Pr\hrd\2 


(22) 


c  = 


&L- 1 


Note  that  c  should  not  be  negative,  i.e.,  c  >  0,  which  implies  that  hrd | 2  <  ^L_A? li)P  ■  Therefore,  the  condi¬ 
tional  probability  (22)  can  be  calculated  as  follows 


QL,k  = 


70 

(L-k-l)Pr 


where 


'|Ndl2=0  J\had\2=0 


Xrd7Q 

Bij  =  1  —  e  — 


Xsde  x*d\hsd\2 Xrde  xrd\hrd\2 d\hsd\2d\hrd\2  A  Bl_ i>fc,  (23) 


Xrd 


Xsd(i-j)Pr  _  X 


Hi 


rd 


For  k  =  L,  the  conditional  probability  c//^  is  given  by 


Arrf70 

e  o-j)pr  —  e 


Xsd7Q 


sd  70 


W,L  =  Pr  [7d,L-i,L-i  <  To]  =  1  -  e  nr-i 


(24) 


(25) 


Based  on  the  above  analysis,  the  average  total  transmission  power  in  (9)  for  the  cooperative  H-ARQ  relaying 
protocol  can  be  given  specifically  in  the  following  theorem. 


Theorem  1  In  the  cooperative  H-ARQ  relaying  protocol,  the  average  total  transmission  power  is 


sd'VO 


L- 1  /l- 1 


t  L—l 


P{L)  =  PsAe  Wr  +  ^  2  Cl’k  +  CP  +  E  Dk  +  Dl 

1=2  \k= 1  /  \fc=l  / 
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(26) 


where 


Ci,k 

Ci, i 


Dk 


Dl 


Ql,kPr  [ Tr  =  k\  Ri>k 

/  _-W 70  _Vr70  ,  ,  _ — . 

(•4|  +A2-A3)[e  '->*  —  e.  mti  )  |  ^  Ps  m  +  (l  -  k)Pr 


,m=  1 


QijPr  [Tr  =  l }  Rij 


sd2 0  ■'sd'TQ 


e  ni  -  e  “i-i  1  -  e 


Asr70 


y  ]  Ps,m > 

m=  1 


QL,kPr  [ Tr  =  fc]  R.L,k 

/  _  AsrTO  _^sr7g 

Sl-ii  (  e  nfc  -  e  nfc-i 


^  ]  Ps,m  +  (L  —  k)Pr 


m=  1 


QL,L,Pr  [Tr  =  L]  Rl,L 


Asr70 
'  Q 


1  _  e  »*-i  1  |  1  -  e  1  ]T  Ps,m, 

m=l 


i  <  A:  <  i  -  i, 


1  <  k  <  L-  1, 


in  which  A\,A2  and  A:>  are  specified  in  (18)-(20),  and  Bl— ] ,k  /.y  specified  in  (24),  respectively.  Also,  70  = 
70A/0,  70  is  the  target  SNR,  A/o  is  the  additive  white  noise  variance,  and  l/\sd,  1/Asr  mid  1/Ard  are  t/i£ 
variances  of  the  source-destination,  the  source-relay  and  the  relay-destination  fading  channels,  respectively.  □ 

When  L  =  2,  the  average  total  transmission  power  expression  P  in  (26)  reduces  to 

P{L  =  2)  =  PS)1  +  (l  -  e  J  {P8, 2  +  e  Pr J .  (27) 


3.1.4  Asymptotically  Tight  Approximation  of  the  Average  Total  Transmission  Power 

The  general  closed-form  expression  in  (26)  for  the  average  total  transmission  power  is  certainly  complicated. 
In  this  section,  we  tty  to  develop  a  simple  and  tight  approximation  of  the  average  total  transmission  power. 

When  L  =  2,  if  we  approximate  e~x  by  1  —  x  for  small  x  (corresponding  to  high  SNR),  then  the  average 
total  transmission  power  P  in  (27)  can  be  approximated  as 


P(L  —  2)  ~  Ps  1  +  (Ps,2  +  Pr 


.Kao 
Ps,  1  ’ 


(28) 


which  is  tight  at  high  SNR.  When  L  >  3,  however,  it  is  difficult  to  derive  the  approximation  of  P(L  >  3) 
directly  from  (26).  So,  we  first  reformulate  the  average  total  transmission  power  P  in  (26)  which  enables  us  to 
develop  a  tight  approximation  at  high  SNR. 

For  simplicity,  let  us  denote  Ci  k  =  7/./,- Pr  [Tr  =  k]  for  k  =  1,  2,  •  •  •  ,1  —  1,  and  Cu  =  qi  i Pr  [Tr  =  l}.  Denote 
Dk  =  qL,kPr  \Tr  =  k]  for  k  =  1,  2,  •  •  •  .  L  —  1,  and  Dl  =  qj.jP'c  [Tr  =  L\.  The  new  notation  coiTesponds  to 

the  notation  Ci,k,  Ci.i,  Dk,  and  Dl  defined  in  Theorem  1.  We  observe  that  C)±  =  Ci.k  ^L=i  Ps,m  +  (l  — 

,  and  Dl  =  DL  J2m=i  ps,m ■  Therefore, 


:,m,  Dk  =  Dk  ^m=l  Ps,m 

the  average  total  transmission  power  in  (26)  can  be  rewritten  by  switching  the  summation  order  (in  terms  of 


k)Pr  ,  Cl, l  =  Cl, l  Em=l  Ps,r 


+  (L  —  k)Pr 
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Ps,m,  1  <  m  <  L,  and  Pr  first)  as  follows 


P(L) 

r  *  _  A  ~  1 

=  P8>1  e  +  + 

;=2  fc=i  fc=i 

L— 1  [L-l  2  L  1  L 

+  yy  -Ps,m  yy  yy  +  yy  +  pS)l  yy 

m=2  L  l=m  k= 1  /c=l  J  fc=l 

"L-l  /-l  L-l 

+Pr  y^y^(t-fc)Q,fc  +  y](L-z)^ 

J=2  fe=l  fc=l 

L 

=  y  ^  Ps,m.Qm  +  PrQr , 

(29) 

where 

Qi  = 

m=l 

e  ps,i  _|_  y^  y^  Qjc  +  yy  z)fc, 

(30) 

Qm  — 

;=2  fe=i  fc=i 

L-l  l  L 

y "  y "  Pi.k  +  y  ]  dl,  2  <  m  <  l — i, 

(31) 

Ql  = 

l=m  k=l  k= 1 

L 

(32) 

Qr  — 

fc=i 

L-l /-I  L-l 

y^y^(/-fc)Q,fe  +  y^(L-z)Lfc. 

(33) 

We  observe  that 

z=2  fc=i  fc=i 

EdTO 

Qi  =  Q2  +  e  Ps’x  , 

(34) 

and  for  any  m  =  2,  3  •  •  •  ,  L  - 

1, 

m 

Q rn  —  Qm+1  ^  ^  ^m,/c  • 

fc=l 

(35) 

In  the  following,  we  calculate  Qrn  for  any  rn  =  1,2 ,■■■  ,L.  Since  Qm  can  be  determined  by  using  the 
previously  determined  Qm+ 1  according  to  (34)  and  (35),  we  derive  Ql,  Ql- i,  ■  ■  •  ,  Q i  in  a  reverse  sequential 
order  as  follows.  First,  Ql  in  (32)  can  be  calculated  specifically  as  follows 


Ql 


L- 1 

r,  L-\,k 

k=  1 


+ 


^sr'TQ 

(Bl~  i,k 


Bl-im  i)  + 


(36) 
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where  B^j  is  specified  in  (24).  For  any  l  =  L  —  1,  L  —  2,  •  •  •  ,  3,  Qi  in  (35)  can  be  calculated  as  follows2 


i- 1 


k=  1 
/-I 


Asr7Q  _^sr7g 


Ql  =  Qi+\  +^2(A!  +  A2-  A3)  e  n*  -  e  n*-i  +  e  n>  -  e  1  -  e  n*-i 


Xsd%  Xsd7  0 


*sr7Q 


=  Q«+i  - 


+ 


=  Qi+ 1  - 


+ 


Lfc=i 

/-l 


Asr7Q  _  Asr7Q 


E/  - — ill  /  _  sa  'u  \  /  - - — ivj 

I  e  -  e  °*-i  +  I  1  -  e  l1-6  °1"1 


Xsd% 


xsr  7Q 


Asr7Q  _  AS7"7Q 


Asr70 
Cl, 


^2  Bi-itk  (  nfc  -  e  )  +  (  1  -  e  n‘-i  )  (  1  —  e  “*-i 

Lfc=l 
l-l 

y~]e  nfc  +  (  1  -  e  ni 

Lfc=l 


Asr7p 


Asd^  0 


1—2 


*sr  7Q 


y~]e  (Bi_ijk  —  Bi_ljk+i)  +  (  1  —  e  n*-1 
Lfc=i 


■srHO 


(37) 


in  which  Ai,A2  and  .4 3  are  specified  in  (1 8)— (20),  and  Bt  j  is  specified  in  (24),  respectively.  When  l  =  L  —  1, 
substituting  Ql  in  (36)  into  (37),  we  have 


Ql-i  =  Ql  — 


+ 


L-  2 


*sr7Q 


e  nk  {BL-  1)fc  —  i?L-i,fc+i)  +  (  1  —  e  !1l-1 

Lfc=l 


AkAq 
ft 


=  0 


L-3 


Asr7Q 


y  e  n*=  {Bl- 2,fc  —  -Bl_2,A;+i)  +  (  1  -  e  ni-2 

Lfc=i 


Asri5o 


L-3 


Asr  70 


Xsd% 
2 


=  e  {BL- 2,fe  -  #L-2,fc+i)  +  (  1  -  e  ni- 
fc=i 

When  l  =  L  —  2,  substituting  the  above  derivation  (38)  into  (37),  we  can  calculate  Ql_2  as 


Ql-2  =  Ql-i  — 


+ 


L-3 


^  e  {Bl- 2,fc  —  77_L_2,fc+i)  +  (  1  -  e  W-2 

Lfc=l 


A.srf7Q 


=  0 


L-4 


*sr7Q 


y  e  {Bl- 3,fc  —  -Bl-3,A;+i)  +  (  1  —  e  n£ 

Lfc=l 


■sd70 


L-4 


Asr70 


y  e  {Bl- 3,k  ~  BL- 3,k+i)  +  (  1  -  e  Ql- 


k= 1 


Asd70 


Assume  that  for  any  l  >  Iq{<  L),  it  is  true  that 

1-2  ,  - 
Asr7Q 


— — ^  ru  /  _ •  sd'YQ  \ 

Ql  =  y  e  (Bi_\  k  -  Bi-ijk+i)  +  (  1  —  e  ni-x  j  ,  l  =  L  —  1,  L  —  2,  •  •  •  ,  Iq,  (40) 


k= 1 


V 


(38) 


(39) 


2To  maintain  the  simplicity  of  the  notation,  we  change  subscript  m  in  (35)  into  l. 
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then  for  l  =  Iq  —  1,  we  have 


Qiq—1  Qlo 


Xsd% 
ft; 


Iq  —  2  , 

^ — _  _  AsrTO 

^  e  (Bi0_itk  -  Bi0_i  k+i)  +  (  1  —  e  nio_1 
Lfc=i 


=  o 


+ 


Iq—S 


^sr  7Q 


yy  e  W  (Bi0—2,k  ~  Bi0_2,k+i)  +  (  1  -  e  SV 
Lfc=l 

Iq—3  \  ~  /  x 

u  AsrTQ  / 


Xsd% 
ft» 


yy  e  (Bi0_2,k  ~  Bi0_2,k+i)  +  (  1  -  e  S2'o_ 

fc=i 


sd'l  0 


(41) 


i.e.,  the  result  in  (40)  is  also  true  for  l  =  Iq  —  1.  Thus,  by  induction  we  can  conclude  that  for  any  l  = 
L  —  1,  L  —  2,  •  •  •  ,  3,  we  have 


1—2 


*sr  7Q 


Qi  =  nk  (Bi_i:k  —  Bi-i^+i)  +  [  1  — e 

fc=i 


(42) 


When  l  =  2,  from  Q3  =  e  ni  (i?2,i  —  -62,2)  +  1  —  e 


.sd'TQ 


we  can  calculate  Q 2  in  the  following 


Q  2  —  Q3  +  C2.1  +  6*2,2 

^sd 


—  Q  3  + 


\  _  \rd^2 

^sd  1\ 


_Asr7o  f  VrfAo  X»di  0  \  /  _  NdlO  _  4sd5 0 

e  ni  (  e  pr  —  e  n2  +  I  e  n2  —  e  ni 


Xsd% 

=  1  —  e  ni  . 


_4m5o 


When  /  =  1,  substituting  Q2  =  1  —  e  ni  into  (34),  we  have 

Asd7Q 

Qi  =  Q2  +  e  =1. 


(43) 


(44) 


We  may  further  approximate  the  probability  Qi  for  any  l  =  2,  3,  •  •  •  ,L.  When  l  =  2,  if  we  approximate  e 
byl  —  2;  for  small  x  (high  SNR),  then  Q2  in  (43)  can  be  approximated  as 


Q2 


-W70 

fli 


(45) 


_  ^sr7 Q 

Note  that  if  we  approximate  e~x  by  1— x+\x2  for  small  x  (high  SNR),  then  for  1  <k<  l— 3,  e 
Si-i^+i)  in  (42)  can  be  approximated  as3 


_ ^sr  7Q 

e  (Bi_i)k  -  Bi_  1^+1) 


^sd'J 0  ^rd'yo  (  1 

20Z_!  Pr  V  -  1  -  k 


1 

l  —  2  — 


(46) 


and  when  k  =  l  —  2, 


_ ^sr7Q 

e  n!-2  [Bi_iti_2  -  Bi_i  i_i) 


^sd'J 0  /  j  ^sr'J 0  .  -VdT 0  \ 

ty-i  v  ty-2  252;.!  +  ) 


(47) 


3When  we  approximate  e  1  by  the  first  three  terms  of  its  Taylor  expansion  1  —  x  +  for  small  x  (corresponding  to  high  SNR), 
then  e~xi  —  e~X2  can  be  tightly  approximated  as  ( xi  —  X2)  [— 1  +  \  (xi  +  *2)]  • 
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Therefore,  for  any  l  =  3, 4,  •  •  •  ,L,Qi  in  (36)  and  (42)  can  be  approximated  as 


Q‘  «  E 


AsrfTO  'VdTO  /  1 _ 1 

2fi{_i  Pr  \l-l-k  l-2-k 


_  Kidl 0  Arrf7)  /  1 _ ,  \ 

2f2z_x  Pr  W-2  v 


Asd7o  /  Asr7o  Artj7o  A  AS(j7o  /  _/W7o  \  A5tj7o  __  1 

^z-i  \^i- 2  2 Pr  )  Qi- 1  l  2fi,_J  J2/-i  2 


/W7o 

Ai/-i 


Asd7o  /  Asr  To  Kd'yo  A 
7v7  2(Z  -  2 )Pr  )  ' 

The  above  calculation  of  the  probability  Qi  can  be  summarized  in  the  following  lemma. 
Lemma  2  The  probability  Qi  in  (30)-(32)  has  a  closed-form  expression  as  follows: 


xsd% 

1  —  e  ni  , 

_  AsrTl 

^1  —  2  „  o. 


X^-=i  e  n,s  (Pz- t,fc  —  S,_i)fc+t)  +  1  —  e  n'-!  I  ,  3  <  l  < 


where  B,  ,  is  specified  in  (24).  Furthermore,  the  probability  Qi  can  be  tightly  approximated  by 


^sd'y 0  (  ^sr'f 0  |  ^rd'y 0 


ffi-1  1  n;_2  ^  2{l-2)Pr  )  ’ 


l  =  2; 


3  <  l  <  L. 


Next,  we  calculate  the  probability  Qr  in  (33).  To  calculate  Qr  more  efficiently,  we  partition  Qr  as  follows 
(corresponding  to  Q2,  Q:\-  •  ■  •  ■  Qi.) 


where 


Qr  —  /*  '  Qr,ii 


Qr, 2  =  Cl,  1  +  D\ , 


L— 2  L-2  /  L— 1 


<3,.  -  EE  Ql.k  +  ^  ]  Pfc  —  E  E  Cz,fc  +  Dk  I  ,  3  <  i  <  L  —  2, 


l=i  k=  1 


k=  1  k=i  \l=k+ 1 


Qr,L—l  ^  ^  C'L— l,fc  H-  ^  ^ 

/c=l  fc=l 

L—l 

Qr,L  =  ^  ^ 


We  observe  that 


Qr,L—l  =  Qr,L  +  ^  C'i-l.fc  —  At,-  1, 
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and  for  any  i  =  3, 4  •  •  •  ,  L  —  2, 


C*  ~  ( 

k=  1  \J: 

In  the  following,  we  calculate  Qr.j  for  any  i  =  2,  3  •  •  •  ,  L.  Since  Q?v,  can  be  determined  by  using  the  previ¬ 
ously  determined  Qr,i+\  according  to  (56)  and  (57),  we  derive  Qv,l,  Qr,L-i,  ■  ■  ■  ,  Qr, 2  in  a  reverse  sequential 
order  as  follows.  First,  Qr,L  in  (55)  can  be  calculated  specifically  as  follows 


i—  1 


L— 1 

E 

=*+i 


Cii  +  77  j 


(57) 


L— 1 

Qr,L  =  ^  '  Bl~ l,k 


k= 1 
L—2 


_  Asr7Q  _  -W7Q 

g  —  g  — 1 


^sr7Q 


*sr70 


As<nO 


=  ]T  e  •«*  (BL_hk  -  Bl_ i)fc+i)  +  e  I  1  -  e  j  , 

fc=i  '  ' 

where  Bi j  is  specified  in  (24).  When  i  =  L  —  1,  Qr,L-i  in  (56)  can  be  elaborated  as  follows 


(58) 


L—2 


Qr,L- 1  =  Qr,L  +  ^  {—Bi,-\,k  +  7?L_2,fc) 
—  Qr,L 


k= 1 
L—2 


xsr  7Q 

flu 


*sr  7Q 


e  nk  -  e  S2fc-i  -  1-e  “i-1  e  “i-1  -  e  W-2 


.2m5o 

T 


xsr 7Q  _  ^sr-7Q 


^sr7Q 


Ksd^  0 


+ 


^  e  nfc  (Bl_ i)fc  -  SL_ljfc+1)  +  (  1  -  e  W-i 
Lfc=l 
L—3 

e  (7?L-2,fc  -  7?L_2,fc+i)  +  e  °L-2  (  1  -  e  n^- 

Lfc=i 


^sr7Q 


AscH0 
2 


(59) 


By  substituting  QT)l  in  (58)  into  (59),  we  have 


L—3 


^sr  70 


Qr,L—i  =  ^  e  (7?L-2,fc  —  7?£,_2,fc+i)  +  e  ni-2  (  1  —  e  nL-2 


fe=i 


-1-e  "i-i 


1-e  W-i 


For  any  i  =  L  —  2,  L  —  3,  ■  ■  ■  ,3,  Qr,i  in  (57)  can  be  calculated  as  follows 

2—1 


Qr,i  —  Qr,i-\- 1  ^  ^  (  ^i,k  “1“  ^i—l,k) 


k= 1 
/  L- 1 


_  Asr7Q  _  ^sr7Q 

g  —  g  — 1 


(~Bii  +  B]_lt)  +  Bl-ij 


_  xsrlQ  _  ^sr7Q 

g  —  g  1 


\/— 2-(-l 
2—1 


^ ^  ,  xsrlQ  xsrl Q  \  f  \srlQ  Xjr 7Q 

Qr,i+i  +  /  (— Bi  k  +  Bi-i  k)  |  e  —  e  nfc_1  j  —  77*^ 

fc=i 


g  _  g  1 


—  Qr,2+1 


+ 


2—1 

x Asr  7Q  AS7-7Q  /  Xsd% 

y^e  nfc  (Bitk  -  Bi;k+i)  +  e  I  1-e  n 
Lfc=i 
2  —  2  , 

x — _  _  Asr 7Q  _ 

e  nfc  (Sj_i)fc  -  Bi- ijfc+i)  +  e  (1-e 

Lfc=l 


scH0 


(60) 


(61) 
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When  i  =  L  —  2,  substituting  Qr,L- 1  in  (60)  into  (61),  we  can  calculate  Qr,L- 2  as 


Qr,L— 2  — 


L— 4 

£< 

fc=l 


*sr  7Q 


^sr7Q 


(-Bl- 3,fe  Bx,-3,A;+l)  +  e  f2L_3  (  1  “  e 


sd^O 


*sr  7Q 


-  1  -  e  "i-i 


A»,j70 

1  -  e  n£-i 


Assume  that  for  any  i  >  io(<  L  —  1),  it  is  true  that 


i—2 


•^sr’70 


Qr,i  —  ^  ^  6 

fc=l 

_>sr2o\  /  Asrf7  0 

-  (  1  -  e  W-i  1(1  —  e  W-1 


AsdTQ 


{Bi-ik  —  -S*_i,fe+i)  +  e  1  (  1  —  e  1 

i  =  L  —  2,  L  —  3, 


*0) 


then  for  z  =  iq  —  1,  we  have 

Qr,io— 1  —  Qr,io  ~ 


+ 


_  A»r 7Q  /  Asd7 0 

e  (-Bi0-1, fe  -  -Bjo-i.fc+t)  +  e  ^o-1  (l— e  ^o-1 

Lfc=l  A 

(£Z?  _  Asr-7n  _  Ast-70  /  _  Asd70 

y~l  e  Qfc  (Bio-2, k  ~  Bio_2,k+i)  +  e  Ql°~2  (  1  —  e  n*o-2 

Lfc=i  A 


io~  3 


Asr7Q 


•^sr70 


E  e  n fc  (Bio-2, k  —  Bi0-2,k+i)  +  e  n*0  2  (l—  e  n‘o 

fc=i 


Ajrf^O 


Asr7Q 


-  1  -  e  W-i  1  -  e  “i-i 


jrfTO 


(62) 


(63) 


(64) 


i.e.,  the  result  in  (63)  is  also  true  for  i  =  io  —  1.  Thus,  by  induction  we  can  conclude  that  for  any  i  = 
L  —  2,  L  —  3,  •  •  ■  ,  3,  we  have 


Qr, 


i 


i-2 

Jfe=l 


Asr7Q 

(Bi-i,k  —  +  e 


1  —  e 


Asr7Q 


1  —  e 


sd^O 


^srTQ 

ni-l 


(65) 


When  i  =  2,  Qr-,2  in  (52)  can  be  calculated  directly  as 


±j — x  \  -  \ 

v ^  _  Asr7Q  _ 2vsr7Q 

Qr,  2  =  /,(— B;;1  +  +  i?L-l,ie  ni 

1=2 

_ Ast-70 

=  -Bi^e  ni 

/  Asd70 

=  (  1  —  e  ni 


(66) 


Furthermore,  we  can  tightly  approximate  the  probability  Qr  t  for  any  i  =  2,  3,  •  •  •  ,  L.  When  i  =  2,  Qr;i  in 
(66)  can  be  approximated  as 


Asd7o  f ,  _  \Sr  To  _  'WytA 

y  fll  2 Hi  J 

^sd'lO 
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(67) 


(68) 


in  which  the  terms  A^70  and  A23^70  in  (67)  are  much  smaller  than  1  at  high  SNR,  which  can  be  ignored.  For  any 
i  =  3, 4,  •  •  •  ,  L  —  1,  in  a  similar  way  as  Qi  in  (48),  Qr^  in  (60)  and  (65)  can  be  approximated  as 


A sd'YO  { ^sr'J 0  Asr70  ^rdf 0  \  Asr70  ^sdf 0 

fij-l  V  ^1-2  fij- 1  2 (i  —  2 )Pr  )  Ql-  1  &L-1 

A sd'JO  {  ^srf 0  ^rdlfO  \ 

\  2(i-2)PTJ  ’ 


(69) 

(70) 


in  which  the  terms  and  As^7°  A^7°  in  (69)  arc  much  smaller  than  and  ^|d7°  ^r7° ,  respectively,  at 
high  SNR,  which  can  be  further  ignored.  When  i  =  L,  (}r.L  in  (58)  can  be  approximated  similarly  as 


s- ,  ^  ^sd'y o  f  \Sr  7o  Ard7o  \ 

^r,L  ~  2(L-2)PrJ  ' 

The  above  calculation  of  the  probability  Qr  can  be  summarized  in  the  following  lemma. 


(71) 


Lemma  3  The  probability  Qr  in  (33)  has  a  closed-form  expression  as  follows: 

L 


Qr  —  ^  '  Qr,ii 


i= 2 


where 


Qr,i  —  \ 


f  •^srT'O 

( 

xsdi0  \ 

e  ni  I 

i  - 

v 

e  nr  j  , 

i  =  2 

■^sr70 

Ast-70  /  Xsdiy 0  \ 

sr^i—2  — 

£fc=i e 

nk 

(Bi~ i,fc  —  +  e  ni~1  (  1  —  e  j 

^sr7Q 


-  1-e  W-i  ill  —  e  sT-i  , 


.id'yo 


(72) 


(73) 


3  <  i  <  L  —  1; 


r  o  ^sr7Q  ^srlQ  /  _  ^scPX). 

Efc=i  e  n,=  {BL-i,k  ~  Bl- i,fc+i)  +  e  (  1  —  e  Q‘L~1  )  ,  i  =  L, 


in  which  Il,  j  is  specified  in  ( 24 ).  Furthermore,  the  probability  Qr  can  be  tightly  approximated  by 

L 


Qr 


A sdlo 


+  E 


AsrfTo  (  AsrTo 


i= 3 


+ 


ArdTO 


(Vi  V  fit- 2  2(i  -  2)P, 


(74) 


Based  on  Lemmas  2  and  3,  the  average  total  transmission  power  P(L)  in  (29)  can  be  approximated  as  follows. 

Theorem  2  In  the  cooperative  H-ARQ  relaying  protocol,  the  average  total  transmission  power  can  be  tightly 
approximated  at  high  SNR  scenario  as 


P(L)  «  PsA  +  (P9> 2  +  Pr) 


Asd7o 


S,  1 


+  ^2  (PS,l  +  Pr 


Asd70 


Asr70 


+ 


Ard70 


2=3 


Em=l  Ps,m  \El-^Ps,m  2(Z  -  2)P, 


(75) 

□ 


3.1.5  Optimal  Transmission  Power  Assignment  for  the  Cooperative  ARQ  Protocol 

In  this  subsection,  we  determine  an  optimal  power  assignment  strategy  for  the  cooperative  H-ARQ  relaying 
protocol  based  on  the  asymptotically  tight  approximation  of  the  average  total  power  consumption  developed  in 
the  previous  section.  We  derive  a  set  of  equations  that  describe  the  optimal  power  level  in  each  (re)transmission 
round  and  enable  its  recursive  calculation  with  fixed  searching  complexity. 
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Before  we  formulate  the  problem  of  finding  the  optimal  transmission  power  sequence,  we  need  to  derive 
the  outage  probability  of  the  cooperative  H-ARQ  relaying  protocol  with  maximum  L  (re)transmission  rounds. 
Let  l>iUl.-j.  denote  the  conditional  probability  that  the  destination  decodes  an  information  packet  unsuccessfully 
after  L  (re)transmission  rounds  given  that  the  event  {Tr  =  k}  occurred.  In  other  words,  is  the  outage 

probability  at  the  destination  in  spite  of  the  fact  that  the  relay  started  forwarding  at  the  (k  +  l)-th  round  for 
any  fc  =  1,2,---  .  L  —  1.  Therefore,  the  outage  probability  of  the  cooperative  H-ARQ  relaying  scheme  after  L 
(re)transmission  rounds  can  be  represented  as 


L 

pout,L  =  J2  p°^k  .  pr  \Tr  =  *.]  7  (76) 

k=\ 

where  the  probability  Pr  [Tr  =  k]  is  specified  in  (12)— (14).  This  outage  probability  is  a  necessary  condition  to 
determine  the  optimal  transmission  power  sequence.  The  conditional  outage  probability  pj!' ’  Lk  in  (76)  can  be 
evaluated  as 

VT*=k  =  Pr  Vld,L,k  <  7o]  >  (77) 

where  the  overall  SNR  7 d,L,k  is  specified  in  (5).  When  Tr  =  L,  we  have 

.  r  _  X.id;y 0 

Prr=L  =  Pr  [7 d,L,L  <  70]  =  1  -  e  .  (78) 

For  any  Tr  =  k,  k  =  1,  2, ...,  L  —  1,  the  conditional  outage  probability  is  given  by 


out.L  ti 

pTr=k  =  Pr 


\hsd\2  < 


70  -  (L  -  k)Pr\hrd\ 

H  T, 


= 


(79) 


where  J  is  specified  in  (24).  Finally,  based  on  the  probability  Pr  [Tr  =  k]  in  (12)— (14)  and  the  conditional 
outage  probability  in  (78)  and  (79),  we  can  obtain  the  outage  probability  for  the  cooperative  H-ARQ 

relaying  scheme  after  L  (re)transmission  rounds  as  follows 


L—l 


P 


out,L  _ 


k= 1 
L—l 


^srl  Q 


y ^BL,k(e  -e  Qfc-1  )  +  (l  —  e 


sd'y  0 


1  -  e 


•^sr7Q 


e  nk  (BL:k  -  BLik+ 1)  +  (  1  -  e  n 

k=\ 


Xsd% 
L 


(80) 


We  can  further  approximate  pout-L  in  (80)  for  high  SNR  scenario  as 


pout,L  ^  ^sd'lo  ^  'WT'O  ^rd'YO 


&L  \^L-1  2(L—l)Pr 


(81) 


In  the  following,  we  determine  an  optimal  power  sequence  P  =  77  j  .  ....  kkjp  P,]  for  the  cooperative  H- 

ARQ  relaying  protocol  such  that  the  average  total  transmission  power  for  the  protocol  to  deliver  an  information 
packet  is  minimized.  We  assume  that  the  relay  transmitted  power  Pr  is  fixed  over  all  retransmission  rounds. 
For  the  H-ARQ  relay  protocol  with  a  targeted  outage  probability  po,  the  problem  of  finding  optimal  power 
assignment  per  round  can  be  formulated  as 


min  P  with  respect  to  Psp .  PS) 2,  •  •  •  ,  Ps,l  ;  Pr  >  0 
subject  to  pout’L  <  p0 


(82) 


where  P  and  poutP  arc  specified  in  (26)  and  (80),  respectively.  Since  solving  the  optimization  problem  with  the 
closed-form  expressions  in  (26)  and  (80)  is  not  analytically  tractable,  we  use  the  asymptotically  tight  approxi¬ 
mation  results  of  the  average  total  transmission  power  and  the  outage  probability,  which  arc  derived  in  (75)  and 
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(81),  respectively.  Then,  the  optimization  problem  is  given  by 


min  P  |  (  p  |  p  \  ^sdTO  \  '  /  p  I  p  \  Md'JO 

;ft>«  Pm +  (p’2  +  p)“or  +  £;  <P'J  +  p)tv7  l  SV7  +  2(T^: 


subject  to 


^sd'J 0  ^sr'J 0 

O.L 


+ 


Ard70 


<  Po 


(83) 


Ml- i  2(L  —  l)Pr 

where  =  Em=i 

Before  we  solve  the  above  optimization  problem,  we  first  show  that  the  average  total  transmission  power 
P  reaches  its  minimum  when  the  constraint  in  (83)  holds  with  equality.  If  there  exists  a  power  sequence 


p*  p* 

/s,l’  ^s, 


2’ 


,  P*L,  P*  such  that 


\sd7o  (  AsrTO 


nt 


+ 


2(L-1)P* 


<  po  where  9*  =  J2m=i  P s,m >  and  the  average 


total  ttansmission  power  P  is  minimized,  then  let  us  consider  another  power  sequence 


Ps,m  —  Ps,m ? 
n.L  =  r/P*L, 
p  —  p* 

1  r  -  J-  r*  i 


1  <  m  <  L  —  1 , 


where  q  is 


A 

V  = 


1 


^sd'YO  /  ^srTO 


n 


p;,l  l  po 

We  can  see  that  the  new  power  sequence  PSii,  Ps2, 

(  ^■sr'YO 


+ 


MdH 0 


-  9 


L—l  2(1/  1)P* 

,  Pr  satisfies 
Md^fO 


L—l 


+ 


V^L-t  2(L  —  1)P, 


=  Po, 


(84) 


(85) 


where  9t  =  X)L=i  Since 


AsdTO  f  'W'TO 


PO 


+ 


^rd^fO 


L—l  2(L  1)P* 

1 


<  so  i]  is  upper  bounded  as 


V  <  ~j5F~  _  ^l-i]  — 


(86) 


s,L 


and  the  corresponding  average  total  transmission  power  is 

P  (Ps, t,  Ps,2i  •  •  •  ,  Ps,L;  Pr)  =  P*1  +  (P*2  +  K) 

L—l 

+  j2(p;,i+p: 


^sdlO 

9* 

.  ^'sd'l 0  ( Asr70 


1=3 


9 


i-i 


9 


+ 


Kdi o 


+  (vPtL+p:)^(^  + 


9 


L- 1 


i— 2  2((  ^)P, 

Ard70 


n 


L—2 


2 (L  -  2)p; 


^  p  /  p*  p*  p*  .  p*\ 

^  S,  2?  5  J-S,L’>  /’ 


(87) 


which  implies  that  the  average  total  transmission  power  resulting  from  the  new  power  sequence  Ps,i,  Ps, 2 ,  •  •  ■  ,  71., 
is  less  than  that  based  on  the  power  sequence  P*l5  P*2,  •  •  •  ,  P*  f .  P*.  This  is  contradictory  to  the  assumption 
that  the  power  sequence  P*  ] .  P*2,  •  •  •  ,  P*L,  P*  minimizes  the  average  total  transmission  power  P.  There¬ 
fore,  the  minimum  average  total  transmission  power  P  can  be  achieved  at  the  boundary  of  the  constraint  (with 
equality)  in  (83). 

Next,  based  on  the  constraint  in  (83)  with  equality,  we  solve  the  optimization  problem  by  considering  a 


Lagrange  multiplier  method.  Denote  F,  =  ,  , 

be  formulated  as 


^sdl' 0  (  I  ^rd'P 0 

fii-l  'r  2(i— l)Pr 


,  then  a  Lagrangian  objective  function  can 


C  (P,  A)  =  PsA  +  (Ps, 2  +  Pr)^  +  ^2(Ps,l  +  Pr)Fl- 1  +  A 


1=3 


Fl  ~Po 


(88) 
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Further  denote  G,  =  Xs^sr^°  and  Mi  =  ^{t-ipp2 »  then  l'lc  derivatives  of  £(P,A)  with  respect  to  PSyT 
1  <  m  <  L,  and  Pr  are 

dC 


dPs,i  ~  1  (^>2  +  Pr) 


^sd'y o 
2 
1 


t  ffi-l+Gl-j  _  \  Fl+Gl 


(  TD  i  ID  \  \  1 

l=3\-^S,l  ~r  K)  A 


ft  =  %?  -  (ft, 3  +  ft)nt  -  EL(P.J  + 

=  Fm_!  -  (P,,m+1  +  -  Ef.m+2(ft,(  +  -  ASt5Ti- 

m  =  3,  4,  •  •  •  ,  L  —  2, 


6>P, 

9£ 


ap, 


dC 


=  FL-2-(Ps,L  +  Pr)^-\ 


,  Pl+Gl 


<9£ 


OP, 


s,L 


=  ft-!  -  Jl 


d£  _  Asd7o 


apr 


r2i 


+  12 1=3  Fl- 1  —  (P«,Z  +  Pr 


fh-i 


—  A^. 

ilL 


Based  on 


d£ 


=  0  and 


dC 


dPs, i  _  w  ap3,2 

<9£  d£ 


=  0,  we  have 


apsa  9ps,2  fii  L 


fil  -  (02  +  Pr)~~^~  -  (Ps,3  +  Pr)Atf"r7°2 


=  0, 


which  implies 


Ps,3  = 


^sd^sr^Q  |_ 


fll  —  (fl2  +  Pr 


'^sd'TO 

fii 


-Pr 


Based  on  -rM—  =  0  and  Af  =  0,  we  have 

Oft,2  oPs,3 


d£  dC 


dPs,2  0Pf 


s,  3 


1 

f^2 


^2^^  -  (fi3  +  Pr)P2  -  (Ps,4  +  Pr)  AsrfAsr7°2n 


fli 


2si3 


=  0, 


which  implies 


Ps,4  = 


2“3 


sd ^  sr  7o  L 


^2^-(^3  +  Pr)P2 


-Pr- 


<9£ 


For  any  m  =  5, 6,  •  •  •  ,  L,  according  to  - - 

dC  dC  1 


=  0  and 


dC 

dPs,m- 1 


=  0,  we  have 


dP.s.m-2  dPs.m-\  “1/n— 2  L 

which  implies 


dm— 2 Pm— 3  ( Urn— 1  "F  Pr)Pm— 2  (Ps,m  "F  Pr) 


AsdAsrT'o 

^m=2^m— 1 


(89) 


(90) 


=  0, 


Prm  — 


dm— 2^m— 1 


9m-2^m—  3  (dm—  1  "F  Pr)Fm— 2 


-Pr. 


(91) 


AsdASr7o 

Moreover,  by  taking  the  derivative  of  £(P,A)  with  respect  to  A  and  setting  it  equal  to  zero,  we  have  the 


constraint  as  Fl  —  po  =  0.  Furthermore,  from  qpl  =  0  and  ,^p  =  0,  we  have  constraints  as  Fl_i  —  A  gF  =  0 

,  ml  _ 

‘  f2z, 


and  A^° 


+  J2h=3  P/- 1  ~  (Ps-,1  +  Fv)q1  ~  —  A^  =  0,  respectively.  Based  on  these  three  constraints,  we 


come  up  with 


^■sd'J 0 

£7i 


£  r 


+  £ 


Z=3  L 


PZ-1  -  (Ps,t  +  Pr 


AFz-t 

fit-t 


FL-iMl 

Po 


(92) 


We  summarize  the  above  discussion  in  the  following  theorem  which  shows  how  to  determine  the  optimal 
transmission  power  sequence  PSji,  Ps,2,  —  ,  Ps,l  and  Pr. 
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Theorem  3  In  the  cooperative  H-ARQ  relaying  protocol,  to  minimize  the  average  total  transmission  power,  the 
optimal  transmission  power  sequence  Ps.\ ,  PSj  2,  •  •  •  ,  PS:l  and  Pr  satisfies  the  following 


and 


where 


Ps,  3  = 


Ps.4  = 


HlH2 


^sd^srf 0  L 
^2^3 


fil  -  (fi2  +  Pr) 


Asrf7o 

Hi 


-P, 


AsdA  sr  7o  L  fii 


^^-(^3  +  ^2 


Psi  = 


H/-2^/-l 

^sd^sr7o 

Asd7Q 

ftl 


^-2-^-3  —  (^i-1  +  Pr)Fl- 2 

Mz-i 


_  p 

t  r? 

—  Pr,  l  =  5, 6,  ■  •  ■  ,  L 


+  £ 


P/-1  —  (-Ps,;  +  Pr 


;=3  L 


0 


;-i 


P  = 


-^sdTo  (  ^srTO 

“ft 


+ 


KdlO 


Mi  = 


FL-iMl 

Po 

^sd^rd'J  0 


(93) 

(94) 

(95) 

(96) 


fii-1  2(*-l)Pr>/’  *  2(i-l)Pr2’ 

in  which  H,;  =  X^m=i  Ps,m»  7o  —  7oM)»  7o  is  ffte  target  SNR,  Afo  is  the  additive  white  noise  variance,  and 
1/A  S(i,  1/Asr  and  1/A  r(j  are  the  variances  of  the  source-destination,  the  source -relay  and  the  relay-destination 
fading  channels,  respectively.  □ 


From  Theorem  3,  we  observe  that  for  any  given  power  Pr,  Fs.i  and  PSj2,  we  can  determine  power  Ps;i-  In 
general,  for  any  l  =  3, 4,  •  •  •  ,L,  with  power  Pr  and  PSj  1,  A5,,  2-  ■  ■  ■ ,  Ps,i- 1-  we  can  determine  power  Ps  [.  Thus, 
for  any  given  initial  power  Pr,  Ps  1  and  Ps  2,  we  can  determine  power  If  j  recursively  for  any  /  =  3.  4,  •  •  •  ,L. 
The  complexity  of  finding  the  optimal  power  sequence  is  reduced  to  that  of  searching  the  initial  power  Ps,i,  Ps>  2 
and  Pr,  i.e.,  searching  over  a  three-variable  space.  We  can  see  that  in  such  a  way,  the  calculation  complexity 
of  finding  the  optimal  power  sequence  for  the  cooperative  H-ARQ  relaying  protocol  is  fixed  regardless  of  the 
maximum  number  of  (re)transmission  rounds  L. 

Especially,  when  L  =  2,  we  can  solve  the  equations  directly  and  obtain  closed-form  results.  In  this  case,  a 
Lagrangian  objective  function  can  be  formulated  as 


P  (P>  A) 


Ps,l  +  (Ps,  2  +  Pr) 


Xsdfo 

Hi 


+  A 


P2  ~Po 


(97) 


The  derivatives  of  C  (P,  A)  with  respect  to  Ps,\,  Ps, 2,  Pr  and  A  arc 


where 


From  !jp  =  0,  we  have 


dC 
d  Ps 


-  —  1  —  (Ps,2  +  Pr 


'“Ef- 


-  A 


F2+G2 
Q2  ? 


dC  _  Nd'ro 


_  \  F2 

AD.r,  > 


dC  _  X sdia  _  \  Mz 
dPr  —  Hi  A  n2  ’ 

dC  _  771 

1 n  ~  ^2  -po, 


K  ax 


F-2  = 


KdlO 

Ho 


( -^srTO  A rdfo  \  ^  Asc;Asr7o  r 

ntr  +  ^r)'  °2  =  — nf- ’  Ma 


^■sd^rd'J 0 

2 P? 


2  Pr2H2 
Ard7oHl 


(98) 
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Based  on  J£-  =  0  and  Jtf  =  0,  we  have 
1  2 


dC  dC  _  \sd^0  pAsdlo  ,G2  _ 

dp^~ dp^~l~ ~^h  ~ {Ps’2+ Pr)^r "  ' 


By  substituting  (98)  into  (99),  we  have 

2A 


-^Pr2  +  PsAPr  +  P,,lPs,2  +  P2S  1  -  =  0. 

Kd  *sdl0 


Based  on  9p  ^  =  0  and  -ffp-  =  0,  we  have 


dC  dC  _  F2-  M2  _ 

—  — A — -  —  0, 


dPs,  2  dPr 


fU 


which  implies 


From  (100)  and  (101),  we  have 


2A 


sr  d2 


'rd 


Pr  +  Ps.lPr  ~~  Ps,l(Ps,l  +  Ps,l)  —  0. 


Ps,  2  =  Ps 


(  Ps,l  A 

,X  \2Asd7o  ) 


p'2 

-Ts,i 

2Asrf7o  ’ 


(99) 


(100) 


(101) 


(102) 


Pi 

where  the  term  s’d ~o  in  (102)  is  much  larger  than  1  at  high  SNR,  in  which  1  can  be  ignored.  By  substituting 
(102)  into  (101),  we  have 


Pr  = 


A  f(j  I’sA 

4Asr 


4AsrPSii 
KdKdi 0 


p. 


S,1 


A/7/ 1 7 1 
A^dA^rTO 


(103) 


where  the  term  \X"\P*^0  in  (103)  is  much  larger  than  1  at  high  SNR,  in  which  all  l’s  can  be  ignored.  By 
substituting  (102)  and  (103)  into  ^  =  0,  we  have 


p3  _ 
Ns,!  ~ 


2AsdA5r70 

po 


1  + 


2A^Asr7o 


1  + 


4Asr-^'s,l 
^sd^rd'y 0 


-  1 


Po 


(104) 


in  which  the  term 


^  4AsrP S|i  ^ 


in  ( 104)  is  much  smaller  than  1  at  high  SNR,  which  can  be  ignored.  Therefore, 


XsdKd% 

the  optimal  transmission  power  PSj  1,  Ps> 2  and  Pr  arc  given  by 


Ps,  1 
Ps,  2 
Pr 


7o 


7o 


7o 


/  2A^Asr  A  3 
\  Po  ) 

(  Kd&Y 

\  2P20  ) 

f ^sd^rd A  2 

V  2p0  y 


(105) 

(106) 
(107) 


.  -1  .  _1  /  1  .  1_  /  1.2 

From  the  above  closed-form  results,  it  is  obvious  that  (p0  )3  <  (p0  U  <  (77  )3  for  a  small  targeted  outage 
probability  pq.  Thus,  we  may  expect  that  Psp  <  Pr  <  Ps< 2.  Under  the  assumption  that  p—  =  j—  =  y-_  =  \ 
and  Afo  =  1,  if  we  set  70  =  10  dB  and  po  =  10-3,  this  expectation  can  be  verified  in  Fig.  5,  which  shows  that 

Ps,  1  <  Pr  <  Ps,2- 


Approved  for  Public  Release;  Distribution  Unlimited. 
23 


Table  1 :  Comparison  of  the  optimal  power  sequence  by  exhaustive  search  (82)  and  closed-form 
results  (105)— (10V)  with  70  =  10  dB,  po  =  10-3,  L  =  2  and  dB  scale. 


Method 

Pa,  1 

Pa,  2 

Pr 

Exhaustive  search  (82) 

21.16 

29.04 

23.42 

(105)-(107) 

21.00 

29.00 

23.50 

3.1.6  Numerical  and  Simulation  Results 

In  this  section,  we  provide  some  numerical  and  simulation  results  to  illustrate  our  theoretical  development  on 
the  optimal  power  assignment  strategy.  In  all  studies,  we  assume  that  the  variances  of  the  channels  hsd,  hsr  and 
hrd  arc  A 1  ^  =  A*  =  =  1  and  the  noise  variance  is  A'o  =  1. 

For  comparison  purposes,  we  derive  the  average  total  transmission  power  for  the  equal  power  assignment 
strategy.  For  a  targeted  outage  probability  po,  the  equal  power  assignment  approach  also  follows  the  constraint 
in  (83)  with  equality  as  follows 


Asd70  (  XSr  To  Xj-d'y 0  \  

VhTi  2(L-i)prJ  ~  Po 


(108) 


Under  the  equal  power  assignment,  the  source  and  relay  transmission  powers  are  fixed  over  all  (re)transmission 
rounds,  i.e.,  =  Pr  =  P  for  1  <  l  <  L.  Thus,  the  equal  power  assignment  is 


P  =  7o 


'■M2 A  sr  +  Kd ) 

2 L(L  -  l)p0 


(109) 


By  substituting  (109)  into  (75),  the  approximation  of  the  corresponding  average  total  transmission  power  is 
given  by 


~2  ^ 

Pequ(L)  «  p  +  2Xsd^  +  ^Yl 


^sd(  2A  sr  +  Kd) 


PU  (Z- l)(Z-2) 


=  70 


'Aad(2A  sr  +  A  rd) 
2 L(L  -  l)p0 


{l  +  2  L(L  —  2)po}  +  2Asd 


(HO) 


In  Figs.  2,  3  and  4,  we  plot  the  average  total  transmission  power  of  the  cooperative  H-ARQ  relaying  protocol 
by  comparing  the  approximation  result  in  (75),  the  exact  closed-form  result  in  (26)  and  the  simulation  result 
for  the  cases  of  L  =  2,  3  and  4,  respectively.  In  this  example,  we  assume  that  the  source  transmission  power 
is  fixed  over  all  (re)transmission  rounds  and  equal  to  the  relay  transmission  power,  i.e.,  Psj  =  Pr  =  P  for 
1  <  l  <  L.  The  targeted  SNR  is  set  at  70  =  10  dB.  Apparently,  we  can  see  that  the  closed-form  result  of  the 
average  total  transmission  power  matches  exactly  with  the  simulation  curve  in  each  figure.  We  also  observe 
that  the  approximation  of  the  average  total  transmission  power  is  loose  at  low  SNR  and  tight  at  high  SNR.  The 
approximation  curve  is  almost  indistinguishable  from  the  exact  closed-form  result  and  the  simulation  curve  for 
SNR  above  15  dB  in  each  figure. 

In  Figs.  5-7,  we  compare  the  optimal  transmission  power  sequence  from  Theorem  3  and  the  exhaustive 
search  result  based  on  the  original  optimization  problem  in  (82)  (without  approximation).  Note  that  in  Fig.  5 
(L  =  2),  we  plot  the  optimal  power  sequence  obtained  from  (105)— ( 107)  instead  of  Theorem  3.  In  these  three 
figures,  we  assumed  that  the  targeted  SNR  is  70  =  10  dB  and  the  required  outage  performance  is  po  =  10-3. 
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Figure  2:  Average  total  power  consumption  per  information  packet  with  Psj  =  Pr  =  P,  1  <  l  < 
2,  70  =  10  dB,  L  =  2. 


Figure  3:  Average  total  power  consumption  per  information  packet  with  Ps  /  =  Pr  =  P,  1  <  l  < 
3,  70  =  10  dB,  L  =  3. 
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Figure  4:  Average  total  power  consumption  per  information  packet  with  Psj  =  Pr  =  P,  1  <  l  < 
4,  70  =  10  dB,  L  =  4. 


The  maximum  number  of  (re)transmission  rounds  is  L  =  2  in  Fig.  5,  L  =  3  in  Fig.  6  and  L  =  4  in  Fig.  7.  We 
can  see  that  the  optimal  transmission  power  values  resulted  from  Theorem  3  (solid  line  with  V)  match  well 
with  the  exhaustive  search  result  from  the  original  optimization  problem  (solid  line  with  ‘o’).  Moreover,  we 
present  the  comparison  of  the  two  solid  lines  numerically  in  Tables  1,  2  and  3  (with  dB  scale).  When  L  =  2, 
Table  1  shows  that  the  optimal  power  sequence  from  the  closed-form  results  (105)-(107)  is  very  close  to  the 
exhaustive  search  result.  For  comparison,  we  also  include  in  the  figures  the  transmission  power  level  of  the 
equal  power  assignment  scheme  (dashed  line  with  ‘o’).  We  observe  that  the  optimal  power  assignment  strategy 
assigns  less  transmission  power  to  the  source  in  the  first  few  (re)transmission  rounds  and  significantly  large 
transmission  power  to  the  source  at  the  last  (L-th)  round  compared  to  the  equal  power  assignment  strategy.  It 
provides  some  insightful  information  on  how  much  power  should  be  assigned  to  the  source  at  each  round  for 
saving  the  average  total  power  cost. 

In  Figs.  8,  9  and  10,  we  plot  the  average  total  transmission  power  for  both  the  optimal  power  assignment 
scheme  and  the  equal  power  assignment  scheme,  with  different  targeted  SNR  70  (from  0  dB  to  25  dB).  The 


Table  2:  Comparison  of  the  optimal  power  sequence  by  exhaustive  search  (82)  and  Theorem  3 
with  70  =  10  dB,  po  =  10”3,  L  =  3  and  dB  scale. 


Method 

PsA 

Ps,2 

Ps,  3 

Pr 

Exhaustive  search  (82) 

16.53 

19.98 

29.65 

17.96 

Theorem  3 

17.20 

20.13 

29.58 

17.63 
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Table  3:  Comparison  of  the  optimal  power  sequence  by  exhaustive  search  (82)  and  Theorem  3 
with  70  =  10  dB,  po  =  10”3,  L  =  4  and  dB  scale. 


Method 

PS,  1 

P.s.2 

Ps,  3 

Ps,4 

Pr 

Exhaustive  search  (82) 

14.39 

14.91 

21.85 

28.44 

15.62 

Theorem  3 

15.54 

15.48 

21.65 

28.30 

15.05 

Transmission  round  index 


Figure  5 :  Transmission  power  sequence  of  the  optimal  and  equal  power  assignment  strategies  with 
70  =  10  dB,  po  =  10"3,  L  =  2. 


required  outage  performance  is  set  at  po  =  10-3.  When  L  =  2,  from  Fig.  8  we  observe  that  the  optimal  power 
assignment  saves  about  1.5  dB  in  average  total  transmission  power  compared  to  the  equal  power  assignment. 
When  L  =  3  and  4,  Figs.  9  and  10  show  that  the  optimal  power  assignment  scheme  outperforms  the  equal 
power  assignment  scheme  with  a  performance  improvement  of  about  2.6  dB.  Moreover,  it  is  interesting  to 
observe  that  in  each  figure,  the  performance  gain  of  the  optimal  power  assignment  scheme  is  almost  constant 
for  different  targeted  SNR  70  (from  0  dB  to  25  dB). 

We  also  compare  the  average  total  transmission  power  required  in  the  two  power  assignment  strategies  with 
different  targeted  outage  probability  values.  We  assume  the  required  SNR  is  70  =  10  dB.  Figs.  11,  12  and  13 
present  comparison  results  for  the  cases  of  L  =  2,  3  and  4,  respectively.  From  the  three  figures,  we  can  see  that 
for  an  outage  performance  of  po  =  10-5,  the  power  savings  of  the  optimal  power  assignment  compared  to  the 
equal  power  assignment  are  5  dB  when  L  =  2,  10  dB  when  L  =  3,  and  1 1  dB  when  L  =  4.  We  also  observe 
that  the  lower  the  required  outage  probability,  the  more  performance  gain  of  the  optimal  power  assignment 
strategy  compared  to  the  equal  power  assignment  strategy. 
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Figure  6:  Transmission  power  sequence  of  the  optimal  and  equal  power  assignment  strategies  with 
70  =  10  dB,  po  =  10"3,  L  =  3. 


Figure  7 :  Transmission  power  sequence  of  the  optimal  and  equal  power  assignment  strategies  with 
70  =  10  dB,  p0  =  lO"3,  L  =  4. 
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Figure  8:  Average  total  transmission  power  of  the  optimal  and  equal  power  assignment  strategies 
withpo  =  10~3,  L  =  2. 


Figure  9:  Average  total  transmission  power  of  the  optimal  and  equal  power  assignment  strategies 
withpo  =  10~3,  L  =  3. 
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Figure  10:  Average  total  transmission  power  of  the  optimal  and  equal  power  assignment  strategies 
withpo  =  10~3,  L  =  4. 


Figure  1 1 :  Average  total  transmission  power  of  the  optimal  and  equal  power  assignment  strategies 
with  70  =  10  dB,  L  =  2. 
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Figure  12:  Average  total  transmission  power  of  the  optimal  and  equal  power  assignment  strategies 
with  70  =  10  dB,  L  =  3. 


Figure  13:  Average  total  transmission  power  of  the  optimal  and  equal  power  assignment  strategies 
with  70  =  10  dB,  L  =  4. 


Approved  for  Public  Release;  Distribution  Unlimited. 
31 


Moreover,  we  show  in  Figs.  8-13  the  average  total  transmission  power  from  exhaustive  search  result  based 
on  the  original  optimization  problem  in  (82).  When  L  =  2,  we  can  see  from  Fig.  8  that  the  average  total 
transmission  power  from  closed-form  results  (105)-(107)  matches  tightly  with  that  from  exhaustive  search 
result.  In  Fig.  11,  we  also  observe  that  the  solid  line  with  V  obtained  from  (105)-(107)  is  tight  below  the 
targeted  outage  probability  of  10-3.  When  L  =  3  (Figs.  9  and  12)  and  L  =  4  (Figs.  10  and  13),  we  can 
see  that  the  average  total  transmission  power  from  Theorem  3  is  close  to  that  from  exhaustive  search  result,  in 
which  the  gap  between  the  two  solid  lines  results  from  the  fact  that  the  solid  line  with  V  is  originated  from  the 
approximation  result  in  (75)  and  the  solid  line  with  ‘o’  is  obtained  from  the  exact  closed-form  result  in  (26). 

3.2  Cooperative  Communication  Protocol  Designs  Based  on  Optimum  Power  and  Time  Allo¬ 
cation 

In  this  subsection,  we  design  and  optimize  cooperative  communication  protocols  by  exploring  all  possible  vari¬ 
ations  in  time  and  power  domains.  First,  we  introduce  a  general  system  model  for  cooperative  communication 
protocol  design  with  flexibility  in  time/power  allocation  and  arbitrary  re-encoding  methods.  Second,  we  study 
the  ideal  cooperative  protocol  and  determine  the  corresponding  optimum  time  and  power  allocations  analyti¬ 
cally.  Third,  we  design  the  practical  cooperative  communication  protocol  based  on  linear  mapping  and  optimize 
the  lineal-  mapping  function  to  minimize  the  outage  probability  of  the  protocol.  Fourth,  extensive  numerical  and 
simulation  studies  are  presented  to  verify  our  theoretical  development.  We  use  the  following  notation  through¬ 
out  this  subsection.  Bold  letters  in  uppercase  and  lowercase  denote  matrices  and  vectors,  respectively.  [■)H, 
det(-)  and  Tr(-)  represent  Hermitian  transpose,  determinant  and  trace  operators,  respectively.  \j  is  an  L  x  L 
identity  matrix,  and  diag(Ai, . . . ,  A/j  is  an  L  x  L  diagonal  matrix  with  diagonal  elements  Ai, . . . ,  A l- 

3.2.1  System  Model 

We  consider  a  cooperative  wireless  network  that  consists  of  one  source,  one  destination  and  one  relay  using  DF 
relaying  protocol.  The  cooperative  strategy  is  described  as  followed  with  two  phases.  In  Phase  I,  the  source 
transmits  information  signal  to  the  destination,  and  the  signal  is  also  received  by  the  relay  as  well.  In  Phase  II,  if 
the  relay  is  able  to  fully  decode  the  information  signal,  it  helps  forwarding  the  information  to  the  destination  via 
certain  re-encoding/transform  methods.  Throughout  this  paper,  we  consider  narrowband  transmissions  in  the 
wireless  network  in  which  channel  between  any  two  nodes  is  subject  to  the  effects  of  frequency  nonselective 
Rayleigh  fading  and  additive  white  Gaussian  noise  (AWGN).  We  assume  that  the  channel  state  information 
(CSI)  is  available  only  at  the  receivers,  not  at  the  transmitters.  Nodes  in  the  network  work  in  a  half-duplex 
mode  where  they  cannot  transmit  and  receive  simultaneously  in  a  same  frequency  band. 

More  specifically,  in  Phase  I,  the  source  broadcasts  its  information  signal  xs  it)  to  both  the  destination  and 
the  relay  during  time  interval  (0,  Tj].  The  received  signals  ys,d(t)  and  ys,r{t )  at  the  destination  and  at  the  relay 
can  be  modeled  as 

ys,d(t)  =  \fP~i  hS:d  xs(t)  +  ns,d(t),  0  <  f  <  Ti,  (111) 

ys,r(t)  =  \fp[  hsr  xs(t)  +  nStr(t),  0  <  t  <  T, ,  (112) 

where  P\  is  the  transmission  power,  xs  (t)  is  normalized  with  power  1,  and  nSjj(t)  and  nsvr  (t)  are  corresponding 
received  noise  at  the  destination  and  the  relay.  In  (1 1 1)  and  (1 12),  and  hs.r  are  the  channel  coefficients  from 
the  source  to  the  destination  and  the  relay,  respectively.  In  Phase  II,  if  the  relay  fully  decodes  the  information 
from  the  source,  then  the  relay  re-encodes  the  information  and  forwards  it  to  the  destination  with  power  P-2, 
otherwise  the  relay  remains  idle.  We  denote  the  time  duration  of  Phase  II  as  T2.  Then,  the  received  signal 
yr,d(t)  at  the  destination  in  Phase  II  can  be  modeled  as 

yr,d(t)  =  \fp-2  K4  xr(t)  +  nr!d(t),  Ti  <  t  <  Ti  +  T2,  (113) 

where  P>  =  Pi  if  the  relay  correctly  decodes  the  information  signal  xs(t),  otherwise  P’i  =  0,  hr(j  is  the  channel 
coefficient  from  the  relay  to  the  destination,  and  nr^{t)  is  received  noise  at  the  destination  in  Phase  II.  In  (1 13), 
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xr(t)  =  A i(xs(t))  is  a  re-encoded  version  of  the  original  information  signal  xs(t)  and  it  is  normalized  with 
average  power  1.  We  note  that  theoretically  arbitrary  re-encoding  function  A4  ( ■ )  may  be  considered  in  the 
protocol  design. 

We  assume  that  the  channel  coefficients  hs^,  hs>r  and  hr(j  are  modeled  as  zero-mean  complex  Gaussian 
random  variables  with  variances  8^d,  52r  and  S2d,  respectively.  The  noise  terms  ns^{t),  ns^r(t)  and  n,. ,/ff) 
are  modeled  as  zero-mean  AWGN  with  variance  Mq.  We  denote  the  total  time  duration  of  each  transmission 
period  as  T  —  Tj  +  T2  and  assume  that  the  fading  channels  are  quasi-static  within  each  transmission  period. 

We  denote  the  ratio  of  the  time  allocation  in  Phase  I  over  the  whole  period  as  a  =  T\/T  €  (0, 1).  If  the  average 
transmission  power  of  the  protocol  in  each  transmission  period  is  P,  then  the  source  and  relay  transmission 
powers  P\  and  P2  should  satisfy  the  constraint 

PiT1+P2T2  =  PT,  (114) 

or  equivalently  aP\  +  (1  —  a)P2  =  P.  For  convenience,  we  further  denote  3  as  the  ratio  of  the  energy  consumed 

in  Phase  I  over  the  total  energy  consumption  in  each  transmission  period,  i.e.,  3  =  so  we  have  3  = 

We  note  that  the  system  model  specified  above  for  cooperative  relaying  protocol  designs  has  two  variable 
domains  to  optimize:  the  power  domain  of  allocating  P\  and  P2  and  the  time  domain  of  allocating  T \  and  T2. 
With  given  power  budget  P,  we  intend  to  optimize  the  power  variables  and  the  time  variables  such  that  the  out¬ 
age  probability  of  the  cooperative  relaying  protocol  is  minimized.  First,  we  study  an  ideal  cooperative  relaying 
protocol  where  there  is  no  constraint  on  the  re-encoding  function  M(-)  and  time  allocation.  The  corresponding 
theoretical  analysis  will  serve  as  guideline  and  benchmark  in  following  practical  cooperative  communication 
protocol  design.  Then,  with  more  realistic  consideration  for  implementation,  we  propose  a  practical  cooperative 
relaying  protocol  by  considering  linear  mapping  technique  as  the  re-encoding  function  M(-),  and  interestingly, 
the  performance  of  the  proposed  practical  protocol  design  is  very  close  to  the  performance  benchmark  of  the 
ideal  cooperative  protocol. 


3.2.2  Ideal  Cooperative  Communication  Protocol 


In  this  section,  we  focus  on  joint  optimization  of  power  allocation  and  time  allocation  in  an  ideal  cooperative 
protocol  where  the  system  can  use  arbitrary  re-encoding  function  A4  ( ■ )  at  the  relay  and  adjust  time  alloca¬ 
tion  arbitrarily  between  Phases  I  and  II.  It  has  been  shown  that  the  relay  and  the  source  can  use  independent 
codebooks  to  achieve  maximum  rate  in  cooperative  communication  protocols  [46].  In  the  ideal  cooperative 
protocol,  we  assume  that  the  source  and  the  relay  perfectly  cooperate  by  using  two  independent  codebooks. 
We  also  assume  that  the  system  can  arbitrarily  allocate  time  duration  to  both  phases,  i.e.,  T\  and  T2  can  be 
arbitrary  positive  numbers  ( T\  +  T2  =  T).  In  the  following,  we  first  calculate  the  outage  probability  of  the  ideal 
cooperative  protocol.  Then,  we  obtain  optimum  power  and  time  allocation  to  minimize  the  outage  probability 
of  the  ideal  cooperative  protocol.  For  any  give  time  allocation  a  €  (0, 1),  we  determine  the  optimum  power 
allocation  at  the  source  and  the  relay  analytically  with  closed-form  expression.  Finally,  we  show  that  in  order 
to  minimize  the  outage  probability  of  the  protocol,  one  should  allocate  more  energy  and  time  to  Phase  I  than 
that  to  Phase  II. 

First,  we  would  like  to  derive  the  outage  probability  of  the  ideal  cooperative  relaying  protocol.  Throughout 
the  paper,  the  outage  probability  is  defined  as  the  probability  that  the  maximum  mutual  information  X  of  the 
transceiver  is  smaller  than  a  predetermined  target  transmission  rate  Ilf  .  If  X  >  Hr,  we  assume  that  the  receiver 
can  decode  the  message  correctly  with  negligible  error  probability. 

In  Phase  I,  with  independently  and  identically  distributed  (i.i.d.)  circularly  symmetric  complex  Gaussian 
input  signals,  the  maximum  mutual  information  between  the  source  and  the  destination  is  given  by 


ls4  =  a  log2 


Pi\hs4\2\ 

A/],  )  ' 


(115) 


in  which  the  time  allocation  ratio  a  shows  the  fact  that  Phase  I  occupies  time  duration  T\  in  each  transmission 
period  (0,  X].  Since  hSJi  ~  CAf(0,  52 d),  then  \hSJi\2  is  an  exponential  random  variable  with  parameter  A.s,;  = 
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1  /Sgd.  Thus,  the  probability  that  the  destination  decodes  incoiTectly  in  Phase  I  can  be  calculated  as 


Pr  [XS)d  <  Rt\  =  1  -  exp  {  - 


M  (n*X. 


2~  -  1 


Similarly,  the  maximum  mutual  information  Is>r  between  the  source  and  the  relay  in  Phase  I  is 


Ps,r  =  a.  log2  1  + 


i  M  |  hs,r 

:  ~ M" 


and  the  outage  probability  that  the  relay  fails  to  decode  the  message  in  Phase  I  is 


Pr  [Xs>r  <  Rt]  =  1  —  exp  <  — 


M0  /  * 


2tt  -  1 


If  the  relay  fully  decodes  the  message  from  the  source,  then  the  relay  forwards  the  message  to  the  destination 
by  using  an  independent  codebook.  With  two  independent  codebooks  at  the  source  and  relay,  the  channels  of 
source-destination  and  relay-destination  can  be  viewed  as  a  pair  of  parallel  channels,  thus  the  joint  maximum 
mutual  information  from  the  source  to  the  destination  in  the  two  phases  is  given  by 


X joint  —  O: 


lo§2  ^  +  (1  -  a)  log2  ^1  + 


Pi  \hr,d\ 

M 


Since  hr>d  ~  CM(0,  S^d),  \ hrJi | 2  is  an  exponential  random  variable  with  parameter  \r  d  =  1  /S^d.  The  proba¬ 
bility  density  function  (pdf)  of  the  random  variable  hr  d 1 2  is  /,  2  (zr.d)  =  M exp  <  —  AT  >  with  zrj  >  0. 

\hr’d\  ’  M  l  °r,d  J 

Similarly,  the  pdf  of  the  random  variable  \hs  d\2  is  given  by  /,  ,2  (zs  d)  =  -J— exp  <  —  >  with  zsd>  0. 

\hs’d\  ’  d*’d  {  j 

Then,  with  the  relay  signal,  the  probability  that  the  destination  decodes  incoiTectly  is  given  by 


r^~ i)  /  r 

Pr  [X joint  <  Rt\  =  J  f\hsdf  (zs,d)  (  J 


,,  ,  iJT-“1°S2(1+plzs,d/A/b)  . 

t(2 - ^ - -1) 


Aq  .  -RT-nl°g2(l+-plz3,dWo) 


exp 


Note  that  the  outage  events  of  the  cooperative  relaying  protocol  have  two  possibilities:  i)  both  the  destination 
and  the  relay  fail  to  decode  the  message  in  phase  I;  and  ii)  the  destination  fails  to  decode  the  message  jointly 
in  Phase  II  even  when  the  relay  fully  decodes  the  message  and  helps  forwarding  the  message  to  the  destination. 
Therefore,  the  overall  outage  probability  Vout  of  the  ideal  cooperative  protocol  is  given  by 


Pout.  =Pr  [XS:d  <  Rt]  Pr  \Xs,r  <  RT\  +  (1  —  Pr  \XS)T  <  RT])  Pr  [Xjoint  <  RT\ 


M  ,n?X. 


M 


„^(2M-l) 


1  —  exp  <  — 


PlS2S}r 


An  (  RT-al°S2(1+plzs,d/'Aro'> 


Mi  * 


=  1  -  exp l  -  7^(2  «  -  1)  \  )  (  1  -  expl  -  -^-(2  j»  -  1)  \]  +  expl  -  7^(2  «  -  1) 


PlS2s,r 
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3.2.3  Optimum  Power  and  Time  Allocations 

In  this  subsection,  we  would  like  to  minimize  the  outage  probability  in  (121)  to  determine  optimum  power  and 
time  allocations  for  the  ideal  cooperative  relaying  protocol.  However,  the  closed-form  expression  in  (121)  is 
not  tractable  for  analytical  purpose.  Therefore,  in  the  following,  we  first  derive  an  approximation  of  the  outage 
probability  (121)  which  is  asymptotically  tight,  then  we  determine  the  optimum  time  and  power  allocations 
based  on  the  asymptotically  tight  approximation  of  the  outage  probability. 

We  would  like  to  use  the  first-order  Taylor  series  approximation,  i.e.  exp(x)  ~  1  +  x  for  x  close  to  0,  to 
simplify  the  two  terms  Pr  [Is,d  <  Rt\  and  Pr  [Ts.r  <  Ilj]  in  (121).  With  high  SNR  where  Pi/Afo  is  large,  we 
have 


Pr  [P-s.d  <  Rt\ 


Pr  [Is.r  <  Rt\ 


W0(2^  -  1) 
P^ld  ’ 
W0(2^  -  1) 


PiSl 


s,r 


Thus,  the  product  of  Pr  [Xs,d  <  Rt]  and  Pr  [ISJ-  <  Rt]  can  be  approximated  as 


Pr  [Ps,d  <  Rt]  ■  Pr  [Xs,r  <  Rt] 


K2 


pfd 


-A(a) 


s,d 


(122) 

(123) 


(124) 


where 


A(a) 


(2»  ~l)s 

S' 2 

us,r 


(125) 


Note  that  A{a)  >  0  and  its  first-order  differential  is  continuous  for  any  a  €  (0, 1).  Furthermore,  we  can  see 
in  (123)  that  the  term  Pr  [XSjr  <  Rt]  «  1  in  the  high  SNR  region,  so  we  have  1  —  Pr  [Is>r  <  Rt]  ~  1  with 
large  Pi  /No. 

In  order  to  obtain  an  asymptotically  tight  approximation  for  the  term  Pr  [ I  joint  <  Rt]  ,  we  need  the  following 
lemma  that  was  developed  in  [40]. 


Lemma  4  [40]  Assume  that  uSl  and  vS2  are  two  independent  scalar  random  variables.  If  their  cumulative 
distribution  functions  (CDF)  satisfy  the  following  properties: 

lim  si  •  Pr[uSl  <  R]  =  a  ■  f(R), 

S 1— »oo 

lim  s 2  •  Pr[vS2  <  R]  =  b  ■  g(R), 

S2 — ^OO 

where  a  and  b  are  constants,  f(R)  and  g{R)  are  monotonically  increasing  functions,  and  the  derivative  of 
g(R),  denoted  as  g'(R),  is  integrable,  then  the  CDF  of  the  sum  of  the  two  independent  random  variables  has 
the  following  property: 

lim  S1S2  • 

51  ,S2 — ^OO 

To  use  Lemma  4,  we  observe  that  Tj)mnt  in  (119)  includes  two  independent  random  variables  which  can  be 
written  as  Tjoint  =  uSl  +  uS2,  where  uSl  =  alog2  ^1  +  and  uS2  =  (1  -  a)  log2 

Let  s  1  =  and  s2  =  and  since  \h.^d\2  and  \hr.d\2  are  exponential  random  variables  with  parameter  1/5% d 


rR 

Pr[uSl  +  vS2  <  R]  =  ab  ■  /  f(r)g'(R  -  r)dr.  (126) 

Jo 
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and  1  /62d,  so  we  have 


lim  si  •  Pr[uSl  <  Rt]  =  lim  si  •  Pr 


S 1— >-oo 


si  — ^oo 


\hs,df  < 


R'-p 

2  2—  -1 


Si 


lim  S2  ■  Pr\uS2  <  Rt]  =  lim  S2  ■  Pr 

S2 — ^OO  S2 — ^OO 


\K,d\  < 


Rrp 

2  2  -  1 


«2 


1  R-t 

TIT  (2—1), 

Ss,d  s - - - ' 

=/(Rt) 


1  . 

=  72- (2^-1) 

°r,d  V ^ - ' 

=9(Rt) 


(127) 


(128) 


We  can  see  that  in  (127)  and  (128),  both  /(Rt)  and  g(Rr)  are  monotonically  increasing  functions,  and  further- 

1  9 

more  we  have  q'(Rt)  =  “  2  i— <*  .  By  applying  Lemma  4,  we  obtain  an  approximation  of  Pr  [I joint  <  Rt]  at 

the  high  SNR  region  as 

fRT 


Pr  [P joint  Rt] 


K2 


PlPoJhS2 


f(r)g'(RT  -  r)dr 


s,ar,d  ^  0 

K2 


'pip^IA 


[Rt  ,  r  ,  In  2  , 

/  (2“  —  1) - 2  ~^Rdr 

Jo  1  —  a 


v„2 


In  2  rt 

2i-a 


(5^(1  -  «)  7o 


fRr 

r(l-2a)  r 

2  o(l-a)  _  2«-l 

Jo 

dr . 


(129) 


=B(a) 


In  (129),  we  observe  that  B(a)  >  0  and  its  first-order  differential  is  continuous  for  any  a  €  (0, 1).  To  further 
simplify  (129),  we  calculate  the  integral  in  B(a)  by  considering  two  cases:  a  =  ]j  and  a  When  a  = 
B(a )  can  be  derived  as 


B(a)  =  2^22Rt  [  T  \  1  -  2-2rl  dr  =  -J-(2RTln2  •  22Rt  -  22Rt  +  1). 
°r,d  «/0  K,d 


(130) 


When  ct/1,  since  fRr  2 crdr  =  ^2  [2cRt  —  l]  for  any  constant  c  >  0,  we  have 


.  In  2  3x_ 

B(a )  =  ,.n  - r21-« 


r Rt  r(l-2a) 

rRr  r 

/  2  “(!-“)  dr  — 

/  2<*~1dr 

Jo 

Jo 

In  2  3t_ 

2  l-a 


6r,Al  ~  «) 


ail  —  a)  (  rt4- 
- 1 - 1 —  I  9  0(1- 

(1  -  2a)  In  2 


-2a) 


2  a(l-a)  —  1  — 


a  —  1  (  3x_ 


In  2 


2^t  -  1 


1 


a  rt  a  —  1  3x_  „  , 

-  2*m  + 1  . 


S2d  U-2 a2  a  +  1-2 a 


(131) 


Therefore,  B(a)  in  (129)  can  be  given  explicitly  as 

f  ^-(2RTln2-22RT -22Rt +  1), 


B(a)  =  { 


a  —  2; 


Rrp 


1_^a  2  «  +  i^2q2  1~Q!  +  1  ,  otherwise. 


(132) 


We  summarize  the  above  discussion  on  the  approximation  of  the  outage  probability  in  the  following  theorem. 


Theorem  4  In  the  ideal  cooperative  protocol,  the  outage  probability  can  be  approximated  as 


Pout 


-  A  MqA(o) 

>  out  — 


X  2  p2 
°s,d 1  1 


,  A r2B(a) 

SldPiPi  ’ 


(133) 


where  A(a)  and  B(a)  are  specified  in  (125)  and  (132),  respectively,  and  the  approximation  is  asymptotically 
tight  at  high  SNR. 
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Figure  14:  Comparison  of  the  exact  and  approximation  of  the  outage  probability  in  three  cases: 
{^rf,  Slr:  8r.d}  =  I1,  h  !}•  {8ldd  8s,r ,  8r,d}  =  U,  10,  and  {8ldd  8s,n  8r.d)  =  {*>  *>  10l-  We  aSSUme  that 

a  =  Rt  =  2,  Mq  =  1,  and  Pi  =  P2  =  P- 

In  Fig.  14,  we  compare  the  exact  value  of  the  outage  probability  which  was  calculated  based  on  (121)  and  the 
asymptotic  approximation  of  the  outage  probability  in  (133)  in  three  scenarios  with  different  channel  variances: 
(i)  {6s,d^Ss,r,Sr,d}  =  {M,l};  (“)  {8‘Lid  6lr :  8r.d}  =  {1,10,1};  and  (iii)  {Sld,6lr,S^d}  =  {1,1,10}-  In 
the  comparison,  we  assume  that  a  =  Rt  =  2,  Mq  =  1,  and  P\  =  R>  =  P.  From  the  figure,  we  can  see 
that  the  approximation  of  the  outage  probability  is  tight  at  reasonable  high  SNR  for  various  channel  conditions. 
The  curve  of  the  outage  probability  approximation  merges  with  the  curve  of  the  exact  calculation  at  an  outage 
probability  of  10-2  in  each  case. 

Next,  we  would  like  to  jointly  optimize  power  and  time  allocations  for  the  ideal  cooperative  protocol  based 
on  the  asymptotically  tight  approximation  of  the  outage  probability  developed  in  Theorem  4.  Note  that  the  time 
allocation  ratio  a  =  T\/T  may  take  any  number  in  the  range  of  (0, 1),  while  the  power  parameters  Pi  and  P2 
should  satisfy  the  power  constraint  in  (114),  i.e.  aPi  +  (1  —  a)  Po  =  P.  The  problem  of  optimizing  time  and 
power  can  be  specified  as  follows: 


^  „  ,  A  NftAlot) 

min  Vout(a,  Pi,  P2)  =  2  2 

0  „ 


Js,d  1 


+ 


A r02B(a) 


a, Pi, P2 

s.t.  aP\  +  (1  —  a)P2  =  P, 

0  <  a  <  1,  Pi  >  0,  P2  >  0. 


8ldP^2 


(134) 


We  found  that  for  any  given  time  allocation  ratio  a  €  (0, 1),  we  are  able  to  express  the  corresponding 
optimum  powers  P\  and  P2  in  terms  of  the  time  allocation  ratio  a  with  closed-form  expressions,  which  are 
denoted  as  P}(a)  and  P2*  (a),  respectively.  Moreover,  we  found  that  for  any  time  allocation  ratio  a  €  (0, 1), 
the  protocol  should  allocate  more  energy  to  Phase  I  than  that  to  Phase  II  (i.e.  /3  =  >  })  in  order  to 

minimize  the  outage  probability  of  the  protocol.  The  results  are  summarized  in  the  following  theorem. 
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Theorem  5  In  the  ideal  cooperative  protocol,  for  any  given  time  allocation  ratio  a  £  (0, 1),  the  corresponding 
optimum  powers  Px*  (a)  and  P|  (a)  are  given  by: 


1  1  +  \Jl  +  8[aA(a)]/[(l  -  a)B(a)]  p 

Oi  3  +  yj\  +  8[a^4(a)]/[(l  —  a)B(a)] 

-1 _  2 

1  —  a  3  +  y7!  +  8[<xA(a:)]/[(l  —  a)B(a)] 


(135) 

(136) 


where  A(a )  and  B(a)  are  specified  in  (125)  and  (132),  respectively.  Moreover,  the  resulting  energy  allocation 
ratio  3  =  2fjApl  [s  strictly  larger  than  which  means  that,  the  protocol  should  allocate  more  energy  to  Phase 
I  than  that  to  Phase  II  in  order  to  minimize  the  outage  probability. 


Proof:  For  simplicity,  in  this  proof  we  drop  the  parameter  a  in  A  (a)  and  11(a),  and  denote  them  as  A 

and  B,  respectively.  Based  on  the  power  constraint  aP\  +  (1  —  a  )  Pi  =  P,  the  relay  transmission  power  If 
can  be  written  as 


P2 


P  —  aP\ 
1  —  a 


(137) 


Since  P>  >  0,  (137)  implies  that  aP\  <  P.  Therefore,  for  any  given  time  allocation  ratio  a  £  (0, 1),  the 
problem  (134)  can  be  reduced  as 


min  n..,i  Pi ) 


MIA  V02(l  -a)B 
SldP?  fldPl(P-aPi) 


s.t.  0  <  a  Pi  <  P. 


(138) 


To  find  the  optimal  power  Pi  in  the  problem  (138),  we  proceed  in  two  steps. 

First,  by  taking  derivative  of  the  target  function  T0ut{Pi )  in  terms  of  Pi,  we  have 

dVout(Pi)  _  A/"02  -2 A(P  -  aP\)2  -  (1  -  a)B(P  -  2aP,)P, 


dPi 


62, 

s,d 


Pf(P-aPi)5 


Let  =  0,  we  have 


— 2y4(P  —  aPi)2  —  (1  —  a)B(P  —  2aP\)Pi  =  0. 

By  solving  the  above  equation,  we  have  two  possible  solutions  as  follows: 

p*  [4aA  —  (1  —  a)B ]  ±  ^/8a(l  —  a)AB  +  (1  —  a)2B2 
1,=l=  4a[aA  —  (1  —  a)B } 


(139) 


(140) 


(141) 


Second,  for  the  two  possible  solutions,  we  would  like  to  show  that  only  the  solution  P* _  can  satisfy  the 
constraint  in  (138).  We  detail  the  discussion  in  two  scenarios: 
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1)  When  aA  >  (1  —  a)B,  in  this  case  the  denominators  of  P‘{  _  and  P*+  arc  positive.  We  have 


Pi 


i,-  • 


P 


i,+ 


aP*  _  = 


> 


> 


[4 aA  -  (1  -  a)B]  -  ^8a(l  -  a)AB  +  (1  -  a)2P2 
4[aA  —  (1  —  a)B] 

[4a  A  —  (1  —  a)B]  —  \J 4a2  A 2  +  4a  (1  —  a)AB  +  (1  —  a)2B 2 
4  [aA  —  (1  —  a)B] 

[4 aA  -  (1  -  a)B]  -  ^J\laA  +  (1  -  a)B]2  _  P 

4[aA  -  (1  -  a)B]  “  ■ 2  ’ 


^  [4aA  -  (1  -  a)B]  -  ^9(1  -  a) 2P2 

^  ir.  1  /i  _  \  nl  J 


aPl+  = 


> 


4[aA  —  (1  —  a)B\ 

[4aA  —  (1  —  a)B]  +  ^/8a(l  —  a)AB  +  (1  —  a)2B2 
4[aA  —  (1  —  a)B\ 

[4 aA  -  (1  -  a)B]  -  ^9(1  -  a)2P2 


P 


4[aA  —  (1  —  a)B\ 


P>  P. 


P 


We  can  see  that  the  solution  P*  +  does  not  satisfy  the  power  constraint  and  only  P*  _  satisfies  the  power 
constraint  in  this  case.  Moreover,  the  corresponding  energy  allocation  ratio  /3  =  "  ^j;  is  within  (^,  1). 

2)  When  aA  <  (1  —  a)B,  in  this  case  the  denominators  of  P*_  and  P*  ,  in  (141)  arc  negative.  For 
convenience,  we  multiply  both  their  denominators  and  numerators  by  —1,  then  we  have 


p* 


P 


i,+ 


«P 


l 


> 

> 

aP*_  < 

aPl+  = 
< 


—  \4aA  —  (1  —  a)B]  +  ^8a{\  —  a)AB  +  (1  —  a)2B2 

—4[aA  —  (1  —  a)B\ 

—  [4aA  —  (1  —  a)B]  +  \J  4a2  A2  +  4a  (1  —  a)AB  +  (1  —  a)2B 2 

—4[aA  —  (1  —  a)B\ 

—  [4aA  —  (1  —  a)B]  +  \J[2aA  +  (1  —  a)B ]2  P 

-4[aA  -  (1  -  a)P]  “  2  : 

—  [4a  A  -  (1  -  a)B]  +  \/ 9(1  —  a)2B2 

— 4[aM  —  (1  —  a)B] 

-[4a A  -  (1  -  a)B]  -  ^8a(l  -  a)AB  +  (1  -  a)2B2 
—4[aA  —  (1  —  a)B] 

—  [4a A  -  (1  -  a)B]  -  (1  -  a)B 

-4[aA  -  (1  -  a)B] 


We  can  see  that  the  solution  Pj*_  is  the  only  feasible  solution  satisfying  the  power  constraint  in  (138), 
and  the  corresponding  energy  allocation  ratio  f3  is  within  the  interval  (5,1). 


The  above  discussion  shows  that  P*  _  is  the  only  solution  that  satisfies  the  power  constraint  and  minimize  the 
outage  probability.  The  resulting  energy  allocation  ratio  (3  is  always  strictly  larger  than  4. 

We  may  rewrite  the  optimum  power  P*  _  as 


p*  _l  +  y/l  +  8(aA)/[(l-a)B]P 
1,_  3  +  yiT8(a]4)7[(r^)Pj  a  ' 

Note  that  A  and  B  arc  shorthands  of  A(a)  and  B(a),  respectively.  For  completeness,  for  any  given  time 
allocation  ratio  a  €  (0, 1),  the  optimum  source  transmission  power  P*  (a)  is  given  by 


^  _  1  +  \J\  +  8[aA(a)]/[{\  —  a)P(a)]  P 

3  +  \J\  +  8[a.A(a)]/[(l  —  a)B(a)]  a 
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(143) 


Based  on  (137),  the  corresponding  optimum  relay  transmission  power  P>  (a)  is  given  by 

2  P 


P2»  = 


(144) 


3  +  iy/l  +  8[ay4(a)]/[(l  —  a)B(ot )]  1  —  a 

Therefore,  we  prove  Theorem  5  completely.  □ 

Based  on  Theorem  5,  we  substitute  the  optimum  power  solutions  P*(a)  and  P2*(a)  'nl°  the  optimization 
problem  (134)  to  find  the  optimum  time  allocation  ratio,  i.e. 

A  MqA(q)  t  MqB(o) 


min  Vout {a)  =  +  .2  o*r  \o*r  \ 

a  dldpt(a)  5s,dp r(a)-P2(«) 

s.t.  0  <  a  <  1. 


(145) 


We  can  see  that  the  optimization  in  (145)  has  only  a  single  variable  a,  and  we  can  apply  numerical  search  of 
the  single  variable  a  over  the  interval  (0, 1)  to  obtain  the  optimum  time  allocation  ratio  a*  that  minimizes  the 
asymptotic  outage  probability  V0ut •  With  the  optimum  time  allocation  ratio  a* ,  based  on  Theorem  5  we  can 
get  the  corresponding  optimum  source  transmission  power  P*(a*)  and  the  optimum  relay  transmission  power 
P2  (a*).  Let  us  denote  (a* ,  P* (a*) ,  P£ (a*))  as  the  solution  of  optimum  power  and  time  allocation  for  the 
ideal  cooperative  protocol.  We  have  the  following  result  regarding  the  optimum  time  allocation  a*. 

Theorem  6  In  the  ideal  cooperative  protocol,  the  optimum  time  allocation  ratio  a*  is  strictly  larger  than 
i.e.  a*  €  i\-  [ ),  which  means  that  the  protocol  should  allocate  more  time  to  Phase  I  than  to  Phase  II  in  order 
to  minimize  the  outage  probability  of  the  protocol. 


Proof:  We  would  like  to  prove  the  result  by  contradiction.  If  there  exists  an  optimum  solution  (a* ,  P* ,  P2 ) 

with  a*  <  })  that  achieves  the  minimum  outage  probability  P*u/  in  the  problem  (134),  we  will  find  another 
solution  that  results  in  smaller  outage  probability  to  contradict  the  assumption. 

With  the  assumption  of  the  optimum  solution  (a*.  P*.  P|)  with  a*  <  the  resulting  minimum  outage 
probability  V*ut  can  be  expressed  as 


V 


* 

out 


SfJP?)2  sbPtP2  ' 


(146) 


~  ^  P*  (l—a*)P*  i 

Let  us  consider  a  new  family  of  resource  allocation  strategy  (q:.  Pi  ,  P2)  =  (ra*,  -f,  y  1  _TaP  )  for  0  <  r  <  W 

We  can  check  that  the  new  resource  allocation  strategy  (d,  Pi ,  P>)  satisfies  the  power  constraint  in  (134). 
Especially  when  r  =  1,  the  new  resource  allocation  solution  is  reduced  to  the  optimum  solution  (a*,P*,P*). 
With  the  new  resource  allocation  strategy,  the  resulting  outage  probability  is 


(A  P  P\  Nq  (rat*)2 A(ra*)  Af02ra* (1  -  Ta*)B(ra*) 
C(r)  -  Vout(a,  Pi,  P2)  -  g2 d{Pta*)2  +  ' 


(147) 


Since  '"2  j(p  2 -  >  dv  —  2  <'  —  1,  we  can  rewrite  the  functions  ,4 (a:)  and  B(a),  defined  in  (125)  and  (129) 
respectively,  as  follows: 


A{a) 
B(a ) 


(In  2)2  I'Rt 

Pir(p  Jo 

(ln2)2 

5r,da(l  ~ 


1 0 


cRt 

2  «  2  “  dvdr 
Rt  pRt~p 


r 

2«  2 1~advdr. 


'o  Jo 


(148) 

(149) 


Thus,  C(t)  can  be  written  as 


C(t )  =  ci 


f *Rt  pRt 


’  T  CX.  /,  TCX 


dvdr  +  C2 


p Rt  pRt~p 


2i  '  dvdr, 
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(150) 


where  Cl  =  and  c2  =  sl^pfapi-a*)  are  Posltlve  constants. 

In  the  following,  we  would  like  to  show  that  the  differential  of  the  function  C{t)  at  r  =  1  is  negative,  i.e., 

9C{r)  | 
dr  I 


AT02(ln2)2 


r_1  <  0.  By  taking  derivative  over  the  function  C'(r),  we  have 


dC(r) 

8t 


=c\  In  2 


(• Rt  pRt 


T=  1 


/  0 


n  4-  r_.  «+r 
- 2 


cr 


+c2a*  In  2 


P-Rt  rRj'—r 


dvdr 


2a *  2'i  [  — 


+ 


(a*)2  (1  —  a*)2 


dvdr. 


(151) 


We  can  see  that  in  (151),  the  first  term  is  strictly  less  than  0.  Next,  we  would  like  to  show  that  the  second 
term  in  (151)  is  non-positive.  We  denote  the  integrand  as  b(v,  r;  a*)  =  2<>*  2  1  ^  ^  j  and  the 

corresponding  integration  domain  as  A  =  {(n,  r)  G  M2.  :  v  +  r  <  Rt},  then  we  can  rewrite  the  second  term 
of  (151)  as  c2cr*  In  2  ffA  b(v,  r;  a*)dvdr.  The  symmetric  property  of  the  integration  domain  A  implies  that  if 
( v ,  r)  €  A,  then  (r,  v)  €  A  as  well.  Thus,  the  second  term  of  (151)  can  be  given  by 


c2cc*  In  2 


b{y ,  r;  a*)dvdr 


c2a*  In  2 
2 


[i b(v ,  r;  a*)  +  6(r,  u;  a*)]dn(ir. 


(152) 


To  prove  the  second  term  of  (151)  is  non-positive,  it  is  sufficient  to  prove  that  the  sum  of  b{v,r\a*)  and 
b(r,v\a*)  is  non-positive  for  any  ( v,r )  G  A.  Since  a*  <  we  can  see  that  only  if  v  <  r,  we  may  have 
b(v,  r;  a*)  >  0,  which  means  that  b(v,  r;  a*)  and  6(r,  v,  a*)  cannot  be  positive  simultaneously.  So,  there  are 
only  two  possible  cases  to  consider.  When  both  b(y ,  r;  a*)  and  b(r,  v ;  a*)  are  non-positive,  the  proof  is  trivial. 
When  either  b(v,  r;  a*)  or  6(r,  v,  a*)  is  positive,  without  loss  of  generality,  we  assume  that  b(v ,  r;  a*)  >  0  and 
b(r,  v\a*)  <  0,  then  we  have 


b(v ,  r;  a*)  +  b(r,  v ;  a*)  =2«*  21-"* 
<2^21^ 


r  +  v 
(a*)2 


r  \ 
(1  -a*)2) 
r  +  v  \ 
(1  —  a*)2  ) 


r  v 

+  2^2!-“* 

<0, 


r  v  \ 

(a*)2  +  (1  -  a*)2) 

(153) 


where  the  first  inequality  is  due  to  the  fact  that  v  <  r.  The  above  result  implies  that  the  second  term  of  (15 1)  is 
non-positive.  Therefore,  we  conclude  that  |t_1  <  0. 

The  result  of  9Cq^'1  |r_1  <  0  implies  that  we  are  able  to  find  a  r  (r  >  1)  such  that  C(r)  <  C(  1).  Since  C(  1) 
is  the  outage  probability  with  the  optimum  allocation  strategy  (a*,  Pj\  P} )  and  C(t)  is  the  outage  probability 
with  another  feasible  allocations  strategy  (ra*,  -  1_rc^*2  ),  the  fact  that  C'(r)  <  C(  1)  for  some  r  >  1 

contradicts  the  assumption  that  there  exists  an  optimum  solution  (a*,  P(.  P})  with  a*  <  Therefore,  we 
prove  Theorem  6  completely.  □ 

Theorem  6  shows  that  the  equal-time  allocation  strategy  in  most  existing  cooperative  relaying  protocol  (allo¬ 
cating  equal  time  among  the  two  phases)  is  not  optimal  in  general. 


3.2.4  Practical  Cooperative  Communication  Protocol  Design  Based  on  Linear  Mapping 

It  is  difficult,  sometimes  may  be  infeasible,  to  implement  the  ideal  cooperative  protocol  as  the  re-encoding 
function  M(-)  and  time  allocation  can  be  arbitrary.  In  this  section,  with  more  realistic  consideration,  we  design 
a  cooperative  communication  protocol  based  on  linear  mapping,  where  the  relay  uses  linear  mapping  as  the 
re-encoding  function  M(-).  Also,  the  practical  cooperative  protocol  allocates  integer  time  slots  in  Phases 
I  and  II.  It  is  much  easier  to  implement  the  linear  mapping  forwarding  method  with  the  time  allocation  of 
integer  time  slots  in  both  phases.  The  theoretical  results  from  the  ideal  cooperative  protocol  in  the  previous 
section  will  serve  as  guideline  and  benchmark  for  the  proposed  linear-mapping  based  cooperative  protocol 
design.  Interestingly,  the  practical  cooperative  protocol  based  on  optimum  linear  mapping  performs  closely  to 
the  performance  benchmark  of  the  ideal  cooperative  protocol.  In  the  following,  we  first  specify  the  practical 
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cooperative  communication  protocol  design  based  on  linear  mapping.  Then  we  optimize  the  linear  mapping  of 
the  protocol  such  that  the  outage  probability  of  the  proposed  protocol  is  minimized. 

We  intend  to  design  a  practical  cooperative  relaying  protocol  with  L  time  slots  in  Phase  I  and  K  time  slots 
in  Phase  II.  We  assume  that  each  time  slot  has  time  duration  Ts.  In  Phase  I,  the  source  broadcasts  a  block  of 
symbols  xs  ( xs  =  (rcs[l], . . . ,  xs[L])T )  using  L  time  slots,  in  which  each  element  of  xs  has  unit  power.  In  Phase 
II,  if  the  relay  successfully  decodes  the  message,  then  the  relay  takes  K  time  slots  to  forward  a  re-encoded 
message  xr  (xr  =  (xy[l], . . . ,  xr[K])T),  where  xr  =  Gx,  and  G  is  a  K  x  L  matrix  representation  of  the  linear 
mapping.  In  the  protocol,  the  total  time  duration  of  Phase  I  is  Tj  =  LTS  and  the  time  duration  of  Phase  II  is 
T2  =  KTS,  thus  the  time  allocation  ratio  a  of  the  protocol  is 


LTS  _  L 
LTS  +  KTS  L  +  K ' 


(154) 


Based  on  the  theoretical  results  in  the  previous  section,  we  know  that  we  should  allocate  more  time  to  Phase  I 
than  that  to  Phase  II  in  order  to  minimize  the  outage  probability  (i.e.,  \  <  a  <  1),  which  means  that  we  should 
choose  K  <  L  in  the  practical  cooperative  protocol  design. 

The  general  system  model  in  Section  3.2. 1  can  be  modified  for  the  linear-mapping  based  cooperative  protocol 
in  a  discrete  version  as  follows: 


y  s,d  = 

y  s,r  = 

\J  P 1  xs  +  n 

\/ P 1  xs  -(-  ns?r, 

(155) 

(156) 

y  r,d  = 

\J P 2  H r,d 

(157) 

where  HSid  =  hS)dJi,  HSjr  =  hs,r^L,  Hr^  =  hr^K,  ys  d  and  y  sr  are  the  signal  vectors  of  size  Lx  1  received 
at  the  destination  and  the  relay  respectively  during  Phase  I,  yr  d  is  the  signal  vector  of  size  K  x  1  received  at 
the  destination  during  Phase  II,  n , and  ns  r  are  the  noise  vectors  of  size  L  x  1  at  the  destination  and  the  relay 
respectively  in  Phase  I,  and  nT)d  is  the  noise  vector  of  size  K  x  1  at  the  destination  in  Phase  II.  In  (155)— (157), 
the  elements  of  the  noise  vectors  n  v/,  ns.r  and  n,.  ,/  are  independent  Gaussian  random  variables  with  mean  zero 
and  variance  A/fj-  The  channel  coefficients  hS)d,  hs,r  and  hrjj  arc  modeled  as  zero-mean  complex  Gaussian 
random  variables  with  variances  8^d,  8gr  and  5^d,  respectively. 

At  the  destination,  it  combines  the  received  signals  ys  d  and  yr  d  from  both  phases  to  jointly  detect  the  original 
message  xs.  The  combined  signal  at  the  destination  can  be  expressed  as 


V  =  f  yS’d  \  =  (  ^Us’dXs  \  A 

d  V  y r,d  J  V  V^H r,dXr  )  V  nr.d  )  ' 

Since  xr  =  Gxs,  the  combined  signal  yd  can  be  further  simplified  as 

y  d  =  +  nd, 


where  H,/  = 


v/^tH  s,d  \  A 
V^Hr,dG  ),nd~ 


ns,d 

H r,d 


and  ~  CAf(0,  A/qIl+A')- 


(158) 


(159) 


3.2.5  Optimum  Linear  Mapping  Design 

In  this  subsection,  we  optimize  the  linear  mapping  matrix  G  of  size  K  x  L  (K  <  L)  such  that  the  outage 
probability  of  the  proposed  linear-mapping  based  cooperative  protocol  is  minimized. 

First,  we  derive  the  outage  probability  of  the  linear-mapping  based  cooperative  protocol.  From  (155),  we 
observe  that  the  channel  from  the  source  to  the  destination  in  Phase  I  is  a  multi-input  multi-output  complex 
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Gaussian  channel.  Thus,  with  i.i.d.  circular  symmetric  complex  Gaussian  inputs,  the  maximum  mutual  infor¬ 
mation  between  the  source  and  the  destination  in  Phase  I  is 


Zs,d  = 


1 


L  +  K 
L 

L  +  K 


log2 


det  f  lL  + 


log2  ^1  +  ^\hs,d\2^j  , 


(160) 


where  the  factor  L  1  K  in  the  first  equality  is  due  to  the  fact  that  the  cooperative  protocol  uses  L  +  K  time  slots 
in  Phase  I.  Similarly,  the  maximum  mutual  information  between  the  source  and  the  relay  in  Phase  I  is  given  by 


T  — 

-L'Sjr  — 


L 

L  +  K 


log2 


1  +  «  r 
■No 


(161) 


In  Phase  II,  if  the  relay  decodes  the  message  from  the  source  correctly  (i.e.,  1+  =  P2),  the  destination  can 
utilize  the  combined  signal  yd  from  two  phases  to  detect  the  original  message  xs.  From  (159),  we  can  see  that 
the  channel  between  the  output  yd  and  the  input  x.s  is  an  equivalent  multi-input,  multi-output  complex  Gaussian 
channel.  Therefore,  the  maximum  mutual  information  of  the  joint  detection  in  Phase  II  is 
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where  the  second  equality  is  due  to  the  fact  that  det  (I l+k  + 


Mo 


) 


det  (IL  +  and  H«Hd  = 


GHG.  Thus,  for  any  given  target  transmission  rate  Rt,  the  outage  probability  of  the 
proposed  linear-mapping  based  cooperative  protocol  is 


Rout  =  Pr  [Ps,d  <  Rt]  Pr  [Ts,r  <  Rt]  +  (1  —  Pr  [Xs,r  <  Rt])  Pr  \Ld  <  Rt]  ,  (163) 


in  which  I S:d,  Ps,r  and  Ld  are  specified  in  (160),  (161)  and  (162),  respectively. 

Next,  we  try  to  optimize  the  linear  mapping  G  in  order  to  minimize  the  outage  probability  in  (163).  We  note 
that  in  (163),  only  the  term  Ld  depends  on  the  linear  mapping  G.  Thus,  to  minimize  the  outage  probability  Vout, 
it  is  equivalent  to  minimize  Pr  [Id  <  Rt]  under  the  power  constraint  Tr  (GGW)  <  K.  The  power  constraint 
ensures  that  the  forwarded  signals  have  unit  average  power.  The  mutual  information  Zd  in  (162)  can  be  further 
calculated  as 
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We  observe  that  'l,i  \ 


does  not  depend  on  the  linear  mapping  G.  Thus,  to  maximize  Xd,  it  is  equivalent  to  maximize  Td„2  under  the 
power  constraint  that  Tr  (GGH)  <  K.  To  further  calculate  Z^2,  we  consider  eigen-decomposition  GGH  = 
U^AU,  where  U  is  a  unitary  matrix,  A  =  diag(Ai, . . . ,  A k)  and  the  eigenvalues  Ai, . . . ,  A k  arc  non-negative. 
Thus,  we  have 
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(165) 


Figure  15:  Outage  probability  of  the  ideal  cooperative  protocol  with  different  time  allocation  ratio  a  €  (0, 1) 
under  the  channel  condition  {5  2  r,  S2d}  =  {1, 1}.  Assume  that  Pi  =  P*  (a)  and  Po  =  P%  (a)  for  any  given 
a  €  (0, 1)  based  on  Theorem  5,  and  P/ No  =  30dB. 

Since  Y^=x  Xi  =  Tr  (GG^)  <  K,  the  problem  of  optimizing  the  linear  mapping  G  is  specified  as 
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We  can  solve  the  problem  by  using  Lagrange  multipliers.  Consider  an  Lagrange  multiplier  //  and  an  Lagrange 
function 
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Differentiating  the  Lagrange  function  with  respect  to  A j,  we  have 
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(167) 


(168) 


(169) 


(170) 


Figure  16:  Outage  probability  of  the  ideal  cooperative  protocol  with  different  time  allocation  ratio  a  £  (0, 1) 
under  the  channel  condition  {<52r,  S^d}  =  {10, 1}.  Assume  that  Pi  =  P* (a)  and  P-2  =  (cr)  for  any  given 
a  £  (0, 1)  based  on  Theorem  5,  and  P/Mq  =  30dB. 


By  substituting  (170)  into  (169),  we  have 


A,;  = 


i  +  PiIM7-M)  +  i  +  Pi\hs,d\2/No 


P2|M7M) 


=  1, 


(171) 


i.e.,  the  optimum  solution  for  the  problem  (166)  is  A*  =  1  for  %  =  1, . . . ,  K.  Therefore,  the  linear  mapping  G 
minimizes  the  outage  probability  in  (163)  if  and  only  if  GGH  =  Ik-  We  summarize  the  above  discussion  in 
the  following  theorem. 

Theorem  7  In  the  linear-mapping  based  cooperative  protocol,  any  K  x  L  (K  <  L)  linear  mapping  G  with 
power  constraint  Tr  ( GG <  K  minimizes  the  outage  probability  of  the  protocol  if  and  only  if  it  satisfies 
GG  h  =  Ik-  Consequently,  we  may  select  any  K  rows  of  an  L  X  L  unitary  matrix  to  form  the  optimum  linear 
mapping  G. 


3.2.6  Numerical  and  Simulation  Results 

In  this  section,  we  provide  some  numerical  studies  and  simulation  results  to  illustrate  the  performance  bench¬ 
mark  of  the  ideal  cooperative  protocol  and  the  performance  of  the  practical  cooperative  protocol  based  on  the 
optimum  linear  mapping.  In  all  numerical  studies  and  simulations,  we  assume  that  the  target  transmission 
rate  is  Rt  =  2  bits/s/Hz,  the  noise  variance  is  A'o  =  1,  and  the  variance  of  the  source -destination  chan¬ 
nel  is  normalized  as  5^d  =  1.  We  consider  three  exemplary  scenarios  with  various  channel  variances:  (i) 

(77  5r,d}  =  {!>  f};  (ii)  {Ss.r;  ^r.d.}  =  I10’  ^  and  (Ui)  {7r=  Sr,d}  =  I1’  10}‘ 

In  Figs.  15-17,  we  plot  the  outage  probability  of  the  ideal  cooperative  protocol  in  terms  of  different  time 
allocation  ratio  a  €  (0, 1)  for  the  three  channel  conditions  (i)-(iii),  respectively.  In  the  study,  we  assume 
P/Mq  =  30dB.  For  any  given  time  allocation  ratio  a  £  (0, 1),  the  corresponding  optimum  power  allocation 


Approved  for  Public  Release;  Distribution  Unlimited. 
45 


Figure  17:  Outage  probability  of  the  ideal  cooperative  protocol  with  different  time  allocation  ratio  a  €  (0, 1) 
under  the  channel  condition  {<52r,  3'2d}  =  {1, 10}.  Assume  that  Pi  =  P* (a)  and  P2  =  /J2* (a)  for  any  given 
a  €  (0, 1)  based  on  Theorem  5,  and  P/Nq  =  30dB. 

P\  =  P* (a)  and  P>  =  P2  (a)  are  calculated  based  on  Theorem  5.  In  each  figure,  we  plot  the  outage  probability 
approximation  based  on  (133)  as  well  as  the  exact  outage  probability  calculated  based  on  (121)  for  comparison. 
In  the  case  of  A2 r.  =  1  and  5^d  =  1,  Fig.  15  shows  that  the  optimum  time  allocation  ratio  is  a*  =  0.66.  The 
corresponding  optimum  power  allocation  is  P*(a*)  =  1.0971P  and  P2* (a* )  =  0.81 15P  based  on  Theorem 
5,  in  which  the  energy  allocation  ratio  3  =  0.7192.  In  case  of  A2r  =  10  and  5^d  =  1,  Fig.  16  shows  that 
the  optimum  time  allocation  ratio  is  a*  =  0.59.  The  corresponding  optimum  power  allocation  is  P}(a*)  = 
1.0196P  and  /J2*  (a*  )  =  0.9718P  based  on  Theorem  5  and  the  energy  allocation  ratio  3  =  0.5969.  In  case  of 
6gr  =  1  and  A2r/  =  10,  Fig.  17  shows  that  the  optimum  time  allocation  ratio  is  a*  =  0.74.  The  corresponding 
optimum  power  allocation  is  P* (a*)  =  1.1112P  and  P|(a*)  =  0.6835P  and  the  energy  allocation  ratio 
3  =  0.8122.  Figs.  15  -  17  show  that  for  all  three  channel  conditions,  the  optimum  time  allocation  ratio  a*  is 
strictly  larger  than  }  which  is  consistent  to  the  theoretical  development  in  Theorem  6,  and  the  corresponding 
energy  allocation  ratio  /3  is  strictly  larger  than  }  which  is  consistent  to  the  result  in  Theorem  5.  Moreover,  we 
observe  that  the  larger  the  ratio  of  the  relay-destination  link  quality  over  the  source-relay  link  quality  (<52d/<52r) 
is,  the  larger  the  optimum  time  allocation  ratio  a*  and  the  corresponding  energy  allocation  ratio  3  are.  In  the 
three  figures,  we  can  see  that  the  approximation  of  the  outage  probability  based  on  (133)  matches  the  exact 
value  of  the  outage  probability  based  on  (121)  very  well. 

Figs.  1 8-20  show  the  outage  probability  performance  of  the  proposed  practical  cooperative  relaying  proto¬ 
col  based  on  the  optimum  linear  mapping  for  three  channel  conditions  (i)-(iii),  respectively.  For  comparison, 
the  figures  also  show  the  performance  of  the  direct  transmissions,  the  performance  of  the  cooperative  protocol 
based  on  equal  time  allocation,  and  the  performance  benchmark  from  the  ideal  cooperative  protocol.  Table  4 
specifies  time  and  power  parameters  (a,  Pi,  )  for  each  cooperative  protocol  under  the  three  channel  condi¬ 
tions,  respectively.  When  ()  2 r  =  1  and  5^d  =  1,  Fig.  18  shows  that  the  proposed  practical  cooperative  protocol 
outperforms  the  equal-time  based  cooperative  protocol  with  performance  improvement  about  1.5dB,  and  the 
difference  between  the  performances  of  the  practical  linear-mapping  based  cooperative  protocol  and  the  ideal 
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Figure  18:  Performance  of  the  proposed  practical  cooperative  protocol  with  the  optimum  linear  mapping  under 
the  channel  condition  {5^  ,  d^d\  =  {1,1}. 


Figure  19:  Performance  of  the  proposed  practical  cooperative  protocol  with  the  optimum  linear  mapping  under 
the  channel  condition  {6^  ,  5^d}  =  {10, 1}. 
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Figure  20:  Performance  of  the  proposed  practical  cooperative  protocol  with  the  optimum  linear  mapping  under 
the  channel  condition  {<52r,  d^d\  =  {1, 10}. 

benchmark  is  less  than  ldB.  When  A2 r  =  10  and  5^d  =  1,  we  can  see  in  Fig.  19  that  there  is  a  2dB  difference 
between  the  proposed  practical  cooperative  protocol  and  the  ideal  performance  benchmark,  but  compared  to 
the  equal-time  based  cooperative  protocol,  the  performance  difference  is  less  than  ldB.  In  this  case,  the  perfor¬ 
mance  gain  over  the  equal-time  based  cooperative  protocol  is  narrowed  because  the  equal-time  allocation  ratio 
is  close  to  the  optimum  time  ratio  which  is  a*  =  0.59  in  this  case.  When  A2r  =  1  and  <52rf  =  10,  Fig.  20  shows 
that  the  performance  of  the  proposed  practical  cooperative  protocol  has  over  2.5dB  improvement  compared  to 
the  equal-time  based  cooperative  protocol.  In  this  case,  the  performance  of  the  proposed  linear-mapping  based 
cooperative  protocol  is  very  close  to  the  ideal  performance  benchmark  (less  than  0.25dB).  From  the  figures, 
we  can  see  that  the  equal-time  based  cooperative  protocol  is  not  optimum  in  general,  and  the  proposed  practi¬ 
cal  cooperative  protocol  with  the  optimum  linear  mapping  has  performance  gain  in  all  three  different  channel 
conditions,  which  varies  according  to  channel  conditions.  Moreover,  we  observe  that  when  the  ratio  of  the 
relay-destination  link  quality  over  the  source-relay  link  quality  (A2d/<52r)  becomes  larger,  the  gap  between 
the  performance  of  the  proposed  practical  linear-mapping  based  cooperative  protocol  and  the  ideal  benchmark 
becomes  smaller. 


Table  4:  TIME  AND  POWER  PARAMETERS  FOR  DIFFERENT  CHANNEL  CONDITIONS 


Kr,^d} 

Resource  Allocation  (a,  Pi,  P2) 

Ideal  Coop. 

Proposed  Coop. 

Equal-time  Coop. 

{1,1} 

(0.66, 1.0971P,  0.8115P) 

(|,  1.0788.P,  0.8423P) 

(1,1.6328P,0.3672P) 

{10,1} 

(0.59, 1.0196P,  0.9718P) 

(|,  0.9948.P,  1.0079P) 

(1,1.2946P,0.7054P) 

{1,10} 

(0.74, 1.1112P,  0.6835P) 

(|,  1.0829P,  0.7513.P) 

(1,1.8560P,0.1440P) 
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4  Results  and  Discussions 


In  this  project,  we  attacked  the  fundamental  problem  of  assigning  optimal  transmission  power  sequence  for 
cooperative  H-ARQ  relaying  protocol  over  quasi-static  Rayleigh  fading  channels.  Finally,  we  arc  able  to  de¬ 
termine  an  optimal  transmission  power  assignment  strategy  for  the  cooperative  H-ARQ  relaying  protocol  over 
quasi-static  Rayleigh  fading  channels  to  minimize  the  average  total  power  consumption.  We  investigated  for 
the  first  time  the  average  total  transmission  power  consumed  by  the  cooperative  H-ARQ  protocol.  Specifically, 
we  developed  an  analytical  approach  to  obtain  a  closed-form  expression  of  the  average  total  transmission  power 
which  is  valid  for  any  maximum  number  of  (re)transmission  rounds  L  allowed  in  the  protocol.  However,  the 
closed-form  expression  of  the  average  total  transmission  power  is  complicated  in  general,  so  we  developed 
a  simple  approximation  of  the  average  total  transmission  power  which  is  tight  at  high  SNR.  Then,  based  on 
the  asymptotically  tight  approximation,  we  determined  the  optimal  power  sequence  that  minimizes  the  aver¬ 
age  total  power  consumption  of  the  protocol  for  any  given  targeted  outage  probability.  We  derived  a  set  of 
equations  that  describe  the  optimal  power  level  in  each  (re)transmission  and  enable  its  recursive  calculation 
with  fixed  searching  complexity.  The  optimal  power  assignment  solution  reveals  that  conventional  equal  power 
assignment  scheme  is  far  from  optimal.  Numerical  and  simulation  results  illustrate  and  validate  our  theoretical 
development.  Specifically,  in  the  first  set  of  numerical  studies,  we  compared  the  average  total  transmission 
power  of  the  cooperative  H-ARQ  relaying  protocol  by  comparing  the  approximation  result  in  (75),  the  exact 
closed-form  result  in  (26)  and  the  simulation  result  for  the  cases  of  L  =  2,  3  and  4,  respectively.  By  compar¬ 
isons,  we  observed  that  the  closed-form  result  of  the  average  total  transmission  power  matches  exactly  with 
the  simulation  curve  in  each  study.  We  also  observed  that  the  approximation  of  the  average  total  transmission 
power  is  loose  at  low  SNR  and  tight  at  high  SNR.  The  approximation  curve  is  almost  indistinguishable  from 
the  exact  closed-form  result  and  the  simulation  curve  for  SNR  above  15  dB  in  each  study.  In  the  second  set  of 
numerical  studies,  we  compared  the  optimal  transmission  power  sequence  from  the  Theorem  and  the  exhaus¬ 
tive  search  result  based  on  the  original  optimization  problem  in  (82)  (without  approximation).  We  observed 
that  the  optimal  transmission  power  values  resulted  from  the  Theorem  match  well  with  the  exhaustive  search 
result  from  the  original  optimization  problem.  An  interesting  observation  is  that  the  optimal  power  assignment 
strategy  assigns  less  transmission  power  to  the  source  in  the  first  few  (re)transmission  rounds  and  significantly 
large  transmission  power  to  the  source  at  the  last  (L-th)  round  compared  to  the  equal  power  assignment  strategy. 
It  provides  some  insightful  information  on  how  much  power  should  be  assigned  to  the  source  at  each  round  for 
saving  the  average  total  power  cost.  In  the  third  set  of  numerical  studies,  we  compared  the  average  total  trans¬ 
mission  power  for  both  the  optimal  power  assignment  scheme  and  the  equal  power  assignment  scheme,  with 
different  targeted  SNR  70  (from  0  dB  to  25  dB).  The  required  outage  performance  is  set  at  po  =  10-3.  When 
L  =  2,  we  observed  that  the  optimal  power  assignment  saves  about  1.5  dB  in  average  total  transmission  power 
compared  to  the  equal  power  assignment.  When  L  =  3  and  4,  we  observed  that  the  optimal  power  assignment 
scheme  outperforms  the  equal  power  assignment  scheme  with  a  performance  improvement  of  about  2.6  dB. 
Moreover,  it  is  interesting  to  observe  that  the  performance  gain  of  the  optimal  power  assignment  scheme  is 
almost  constant  for  different  targeted  SNR  70  (from  0  dB  to  25  dB).  In  the  fourth  set  of  numerical  studies,  we 
compared  the  average  total  transmission  power  required  in  the  two  power  assignment  strategies  with  different 
targeted  outage  probability  values  for  the  cases  of  L  =  2,  3  and  4,  respectively.  From  the  comparisons,  we 
observed  that  for  an  outage  performance  of  po  =  10  A  the  power  savings  of  the  optimal  power  assignment 
compared  to  the  equal  power  assignment  are  5  dB  when  L  =  2,  10  dB  when  L  =  3,  and  1 1  dB  when  L  =  4.  We 
also  observed  that  the  lower  the  required  outage  probability,  the  more  performance  gain  of  the  optimal  power 
assignment  strategy  compared  to  the  equal  power  assignment  strategy. 

While  most  existing  works  on  cooperative  relaying  protocol  designs  considered  equal-time  allocation  sce¬ 
nario,  i.e.  equal  time  duration  is  assigned  to  each  source  and  each  relay,  in  this  project  we  are  able  to  design  and 
optimize  cooperative  communication  protocols  by  exploring  all  possible  variations  in  time  and  power  domains. 
First,  we  analyzed  and  optimized  the  ideal  cooperative  communication  protocol  where  the  system  can  use  arbi¬ 
trary  re-encoding  function  at  the  relay  and  adjust  time  allocation  arbitrarily  between  Phases  I  and  II.  Based  on 
the  asymptotically  tight  approximation  of  the  outage  probability,  we  obtained  the  optimum  strategy  of  power 
and  time  allocations  to  minimize  the  outage  probability  of  the  ideal  cooperative  protocol.  For  any  given  time 
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allocation  a  €  (0, 1),  we  determined  the  corresponding  optimum  power  allocation  at  the  source  and  the  relay 
analytically  with  a  closed-form  expression.  We  also  showed  theoretically  that  in  order  to  minimize  the  outage 
probability  of  the  protocol,  one  should  always  allocate  more  energy  and  time  to  Phase  I  than  that  to  Phase  II  in 
the  protocol.  We  note  that  in  the  ideal  cooperative  protocol  in  which  there  is  no  constraint  on  the  re-encoding 
methods  and  time  allocation,  it  may  not  be  easy/feasible  to  implement  it  in  practical  systems.  Therefore,  with 
more  realistic  consideration,  we  proposed  a  practical  cooperative  relaying  protocol  design  based  on  linear  map¬ 
ping,  where  the  protocol  considers  linear  mapping  forwarding  method  at  the  relay  and  uses  integer  time  slots 
in  Phases  I  and  II.  The  theoretical  results  from  the  ideal  cooperative  protocol  served  as  guideline  and  bench¬ 
mark  in  the  practical  cooperative  protocol  design.  We  also  developed  an  optimum  linear  mapping  to  minimize 
the  outage  probability  of  the  linear-mapping  based  cooperative  protocol.  Extensive  numerical  and  simulation 
results  show  that  the  practical  cooperative  relaying  protocol  based  on  the  optimum  linear  mapping  outperforms 
the  existing  cooperative  protocol  with  equal  time  allocation,  and  more  interestingly,  the  performance  of  the 
practical  linear-mapping  based  cooperative  relaying  protocol  is  close  to  the  performance  benchmark  of  the 
ideal  cooperative  protocol.  More  specifically,  in  numerical  studies  and  simulations,  we  consider  three  exem¬ 
plary  scenarios  with  various  channel  variances:  (i)  {A2,,,  6^d\  =  {1, 1};  (ii)  {()2,,,  dfd}  =  {10, 1};  and  (iii) 
{A2r,  5^d}  =  {1, 10}.  The  variance  of  the  source -destination  channel  is  normalized  as  d’~(/  =  1.  In  one  set  of 
numerical  studies,  we  compared  the  outage  probability  of  the  ideal  cooperative  protocol  in  terms  of  different 
time  allocation  ratio  a  €  (0, 1)  for  the  three  channel  conditions  (i)-(iii),  respectively.  In  the  study,  we  compared 
the  outage  probability  approximation  based  on  (133)  as  well  as  the  exact  outage  probability  calculated  based  on 
(121).  In  the  case  of  <52r  =  1  and  5^d  =  1,  it  shows  that  the  optimum  time  allocation  ratio  is  a*  =  0.66.  The 
corresponding  optimum  power  allocation  is  P*  (a*)  =  1.0971P  and  P|(a*)  =  0.8115P  based  on  Theorem  5, 
in  which  the  energy  allocation  ratio  3  =  0.7192.  In  case  of  A2 r  =  10  and  S^d  =  1,  it  shows  that  the  optimum 
time  allocation  ratio  is  a*  =  0.59.  The  corresponding  optimum  power  allocation  is  P(  ( a* )  =  1.0196P  and 
P*(a*)  =  0.9718P  based  on  Theorem  5  and  the  energy  allocation  ratio  3  =  0.5969.  In  case  of  ci2r  =  1  and 
5 =  10,  it  shows  that  the  optimum  time  allocation  ratio  is  a*  =  0.74.  The  corresponding  optimum  power 
allocation  is  P^(a*)  =  1.1112P  and  P|(cr*)  =  0.6835P  and  the  energy  allocation  ratio  f3  =  0.8122.  We  ob¬ 
served  that  for  all  three  channel  conditions,  the  optimum  time  allocation  ratio  a*  is  strictly  larger  than  }  which 
is  consistent  to  the  theoretical  development  in  Theorem  6,  and  the  corresponding  energy  allocation  ratio  3  is 
strictly  larger  than  {  which  is  consistent  to  the  result  in  Theorem  5.  Moreover,  we  observed  that  the  larger  the 
ratio  of  the  relay-destination  link  quality  over  the  source-relay  link  quality  {5^d/5^r)  is,  the  larger  the  optimum 
time  allocation  ratio  a*  and  the  corresponding  energy  allocation  ratio  (3  are.  In  another  set  of  numerical  studies, 
we  observed  the  outage  probability  performance  of  the  proposed  practical  cooperative  relaying  protocol  based 
on  the  optimum  linear  mapping  for  three  channel  conditions  (i)-(iii),  respectively.  For  comparison,  we  also 
simulated  the  performance  of  the  direct  transmissions,  the  performance  of  the  cooperative  protocol  based  on 
equal  time  allocation,  and  the  performance  benchmark  from  the  ideal  cooperative  protocol.  In  comparisons,  we 
observed  that  when  i  2r  =  1  and  6^d  =  1,  the  proposed  practical  cooperative  protocol  outperforms  the  equal¬ 
time  based  cooperative  protocol  with  performance  improvement  about  1.5dB,  and  the  difference  between  the 
performances  of  the  practical  linear-mapping  based  cooperative  protocol  and  the  ideal  benchmark  is  less  than 
ldB.  When  ci2r  =  10  and  5^d  =  1,  we  observed  that  there  is  a  2dB  difference  between  the  proposed  practical 
cooperative  protocol  and  the  ideal  performance  benchmark,  but  compared  to  the  equal-time  based  cooperative 
protocol,  the  performance  difference  is  less  than  ldB.  In  this  case,  the  performance  gain  over  the  equal-time 
based  cooperative  protocol  is  narrowed  because  the  equal-time  allocation  ratio  is  close  to  the  optimum  time 
ratio  which  is  a*  =  0.59  in  this  case.  When  r)  2r  =  1  and  h2(/  =  10,  we  observed  that  the  performance  of  the 
proposed  practical  cooperative  protocol  has  over  2.5dB  improvement  compared  to  the  equal-time  based  coop¬ 
erative  protocol.  In  this  case,  the  performance  of  the  proposed  linear-mapping  based  cooperative  protocol  is 
very  close  to  the  ideal  performance  benchmark  (less  than  0.25dB).  From  the  comparisons,  we  can  conclude  that 
the  equal-time  based  cooperative  protocol  is  not  optimum  in  general,  and  the  proposed  practical  cooperative 
protocol  with  the  optimum  linear  mapping  has  performance  gain  in  all  three  different  channel  conditions,  which 
varies  according  to  channel  conditions. 
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5  Conclusions 


In  this  report,  we  summarized  our  findings  resulting  from  the  project  and  described  in  details  the  system  models 
and  the  proposed  methods  and  procedures.  First,  we  determined  an  optimal  transmission  power  assignment 
strategy  for  the  cooperative  H-ARQ  relaying  protocol  over  quasi-static  Rayleigh  fading  channels  to  minimize 
the  average  total  power  consumption.  We  investigated  for  the  first  time  the  average  total  transmission  power 
consumed  by  the  cooperative  H-ARQ  protocol.  Specifically,  we  developed  an  analytical  approach  to  obtain 
a  closed-form  expression  of  the  average  total  transmission  power  which  is  valid  for  any  maximum  number  of 
(re)transmission  rounds  L  allowed  in  the  protocol.  However,  the  closed-form  expression  of  the  average  total 
transmission  power  is  complicated  in  general,  so  we  developed  a  simple  approximation  of  the  average  total 
transmission  power  which  is  tight  at  high  SNR.  Then,  based  on  the  asymptotically  tight  approximation,  we 
determined  the  optimal  power  sequence  that  minimizes  the  average  total  power  consumption  of  the  protocol 
for  any  given  targeted  outage  probability.  We  derived  a  set  of  equations  that  describe  the  optimal  power  level 
in  each  (re)transmission  and  enable  its  recursive  calculation  with  fixed  searching  complexity.  The  optimal 
power  assignment  solution  reveals  that  conventional  equal  power  assignment  scheme  is  far  from  optimal.  For 
example,  for  a  targeted  outage  performance  of  po  =  10~3  and  maximum  number  of  (re)transmissions  L  =  4, 
the  optimal  power  assignment  scheme  saves  about  2.6  dB  in  the  average  total  power  consumption  compared  to 
the  equal  power  assignment  scheme.  We  also  observe  that  the  lower  the  required  outage  probability,  the  more 
performance  gain  of  the  optimal  power  assignment  scheme  comparing  to  the  equal  power  scheme. 

Furthermore,  we  designed  and  optimized  cooperative  relaying  protocols  by  exploring  possible  variations  in 
time  and  power  domains.  First,  we  analyzed  and  optimized  the  ideal  cooperative  communication  protocol 
where  the  system  can  use  arbitrary  re-encoding  function  Ad  (•)  at  the  relay  and  adjust  time  allocation  arbitrarily 
between  Phases  I  and  II.  Based  on  the  asymptotically  tight  approximation  of  the  outage  probability,  we  obtained 
the  optimum  strategy  of  power  and  time  allocations  to  minimize  the  outage  probability  of  the  ideal  cooperative 
protocol.  For  any  given  time  allocation  a  €  (0, 1),  we  determined  the  corresponding  optimum  power  allocation 
at  the  source  and  the  relay  analytically  with  a  closed-form  expression.  We  also  showed  theoretically  that  in 
order  to  minimize  the  outage  probability  of  the  protocol,  one  should  always  allocate  more  energy  and  time  to 
Phase  I  than  that  to  Phase  II  in  the  protocol.  We  note  that  in  the  ideal  cooperative  protocol  in  which  there  is 
no  constraint  on  the  re-encoding  methods  and  time  allocation,  it  may  not  be  easy/feasible  to  implement  it  in 
practical  systems.  Therefore,  with  more  realistic  consideration,  we  proposed  a  practical  cooperative  relaying 
protocol  design  based  on  linear  mapping,  where  the  protocol  considers  linear  mapping  forwarding  method  at  the 
relay  and  uses  integer  time  slots  in  Phases  I  and  II.  The  theoretical  results  from  the  ideal  cooperative  protocol 
served  as  guideline  and  benchmark  in  the  practical  cooperative  protocol  design.  We  also  developed  an  optimum 
lineal-  mapping  to  minimize  the  outage  probability  of  the  linear-mapping  based  cooperative  protocol.  Simulation 
results  show  that  the  practical  cooperative  relaying  protocol  based  on  the  optimum  linear  mapping  outperforms 
the  existing  cooperative  protocol  with  equal  time  allocation,  and  more  interestingly,  the  performance  of  the 
practical  linear-mapping  based  cooperative  relaying  protocol  is  close  to  the  performance  benchmark  of  the  ideal 
cooperative  protocol.  We  observed  that  when  the  ratio  of  the  relay-destination  link  quality  over  the  source-relay 
link  quality  (6^ d/ 6 |r)  becomes  larger,  the  gap  between  the  performance  of  the  proposed  linear-mapping  based 
cooperative  protocol  and  the  ideal  benchmark  becomes  smaller. 
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ACRONYMS 


ACK 

ACKnowledgement 

AF 

Amplify-and-Forward 

ARQ 

Automatic-Repeat-reQuest 

AWGN 

Additive  White  Gaussian  Noise 

CDF 

Cumulative  Distribution  Function 

CSI 

Channel  State  Information 

DF 

Decode-and-Forward 

H-ARQ 

Hybrid  Automatic-Repeat-reQuest 

LTE 

Fong  Term  Evolution  (FTE)-Advanced  standard 

MIMO 

Multiple-Input-Multiple-Output 

MRC 

Maximal-Ratio-Combining 

NACK 

Negative-ACKnowledgement 

QoS 

Quality-of-Service 

SER 

symbol  error  rate 

SNR 

Signal-to-Noise  Ratio 
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